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SUMMARY 


A simulation study was conducted to evaluate the delivery performance of an 
earlier proposed fixed-path metering and spacing (M & S) system and to determine 
if use of the microwave landing system (MLS) improved the system’s delivery per- 
formance over that achieved with the current instrument landing system (ILS) . 

The M & S algorithms used in the study are described with additional detail and 
flew charts included in the appendixes. The M & S system assumes an aircraft 
two-dimensional area navigation (RNAV) capability, and the three navigation sys- 
tems compared were as follows: (1) Very high frequency omnidirectional range/ 

distance measuring equipment (VOR/DME) and ILS; (2) VOR/DME and ±40^ MLS; and 
(3) VOR/DME and ±60^ MLS. The study was made with a computerized Terminal Area 
Air Traffic Model (TAATM), which contains pertinent en route effects, and the 
results of the study were compared to a full-crew, pilot- in- the- loop cockpit 
simulation. 

The batch (fast-time) simulation results show that the aircraft 
interarrival-error standard deviation at the outer marker is only reduced from 
9.2 to 7.9 sec when the ±60^ MLS replaces the VOR/DME for navigation during 
final-spacing path-adjustment maneuvers. Though true for the particular geome- 
try studied, further investigation showed this result reflects the interaction 
of two types of approach route. For only straight-in routes there was no veri- 
fied difference between the inter arrival-error standard deviation of the MLS 
and ILS systems; however, there was a 2.6-sec advantage for the MLS system with 
a geometry containing only direct-course-error (DICE) turn final- approach 
routes. For the system studied the wider coverage ±60^ MLS offers no advantage 
over the basic ±40^ MLS system. 

Real-time cockpit data indicate that moving the gate closer to the runway 
threshold, although theoretically possible because of MLS precision, may affect 
the interarrival delivery precision of the M & S system. Interarrival delivery- 
error deviations achieved during the study were substantially less than those 
usually quoted for current manual vectoring procedures. The study also demon- 
strated that most of the terminal area holding is transferred to more fuel- 
efficient, higher altitude en route delay when the en route time-base metering 
is coupled to the terminal scheduling and spacing process in the M & S system. 

A reduction in terminal delay not only saves fuel but increases safety by reduc- 
ing congestion and lowering the work load of both pilots and controllers. 


1.0 INTRODUCTION 

A study was undertaken in response to a Federal Aviation Administration 
(FAA) need to more completely evaluate a particular early-design fixed- path 
metering and spacing (M & S) system (refs. 1 to 4) and to determine if use of 
the microwave landing system (MLS) for approach guidance improved the system’s 
delivery performance. The proposed fixed-path M & S system emphasizes the use 
of RNAV paths together with computer-derived indicated-airspeed commands to the 
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pilot in order to minimize route variations normally used in achieving a precise 
scheduled landing time- The regulation of the flow rate (metering) into the 
terminal is accomplished via a time-based process - 

Limited simulation evaluations (refs. 1 to 4) of a metering and spacing 
system using speed control, transition vectoring, and two-dimensional RNAV with 
direct engage capability were encouraging. These studies demonstrated the work- 
ability of this type of system and yielded results comparable to M & S systems 
using only vectoring. Another evaluation, which used a single-aircraft root- 
sum-square error analysis (ref. 5) , calculated the delivery performance of a 
strictly speed-control (fixed-path only) system and also that of a speed-control 
system with fine tuning through RNAV direct-engage time-to-turn maneuvers in the 
base-leg region. This analysis indicated that a wide-coverage MLS would have 
significant superiority over other navaids in its ability to support a fixed- 
path speed-control system. 

A more rigorous investigation through simulation was deemed necessary to 
consider the system performance of the MLS fixed-path system. A more realistic 
and detailed simulation would evaluate the multiple aircraft interactions of 
the system subjected to various traffic demands. Also needed was verification 
of theoretical results on such man-machine issues as communication delays and 
pilot reaction under realistic work-load conditions. 

Facilities used in the investigation are the Terminal Area Air Traffic 
Model (TAATM) and the Terminal Configured Vehicle (TCV) Aft Flight Deck Simu- 
lator at the Langley Research Center (ref. 6). The fixed-path M & S control 
logic was incorporated into TAATM and its geometry configured to the Denver 
Terminal Area used in references 1 to 5. Current operational geometry in Denver 
is different, but in order to make comparisons with previous results, the origi- 
nal configuration was used in this investigation. The TAATM is a flexible 
dynamic model of the airborne, navaid, ground control, and communications 
aspects of the terminal area environment which can run in either fast time 
(batch mode) or in real time. The TAATM also models pertinent en route effects. 
Pilot-in- the-loop experiments were conducted with the TAATM and the TCV simulator 
integrated together. The TAATM in conjunction with a simulation operator pro- 
vided an operational environment and voice commands to the TCV simulator. The 
TCV simulator then became one of the aircraft to be controlled in the TAATM 
airspace. 

The report is organized so that the system studied is first described and 
then the investigation results and conclusions follow. Sections 2.0 and 3.0 
explain the terminal geometry and the fixed-path, speed-control M & S system 
evaluated. Section 4.0 describes the models used to represent navigation and 
system errors. Limitations in the M & S system were observed during the inves- 
tigation and presented in section 5.0. Section 6.0 presents the fast- time 
(batch) study results whereas section 7.0 presents the real-time pilot-in-the- 
loop findings. The entire investigation results are summarized in section 8.0. 
Appendix A gives the derivation of a procedure for determining system inter- 
arrival time-error performance from single-aircraft, route-dependent arrival 
errors. This allows a direct comparison of real-time results with those 
obtained from system inter arrival-error performance measured in fast- time runs. 
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/^pendixes B and C give more details of the M & S system and aircraft navigation 
modes and appendix D contains M & S control-logic flow diagrams. 

Several of the figures are exact duplicates of the displays used, and as 
such they contain measurements in feet. To convert these values to meters, use 
the following conversion: 1 ft = 0.3048 m. Also, knots are indicated in sane 

figures as K. 

Use of trade names or names of manufacturers in this report does not con- 
stitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


2.0 AIR TRAFFIC CONTROL (ATC) TERMINAL AREA METERING 
AND SPACING ENVIRONMENT FOR DENVER 

The terminal area environment for the metering and spacing studies was 
designed based on the proposed specifications in reference 4. The overall con- 
figuration assuming MLS capability is depicted in figure 1 , which shows the 
merger of two of the seven arrival routes into a single straight-in approach 
from the east (BYERS), merger of two routes into a north corner-post approach 
(LONQIONT) , one southwest corner-post approach (SHAWNEE) , and merger of the 
last two routes into a south corner-post approach path (ELIZABETH) . The four 
ATC intersections (VOR/DME defined) for which the standard terminal arrival 
routes (star's) are named are the so-called metering or initial arrival fixes 
(lAF's) at which holding stacks are deployed when necessary. Each path fran 
the terminal area perimeter to the final-approach gate is a two-dimensional area 
navigation (RNAV) path defined by navigation waypoints. Exceptions to this 
are the racetrack-shaped holding patterns and the delay vector areas. The way- 
point names and locations do not necessarily represent reality in the Denver 
area but were chosen to satisfy the experimental requirements for the TAATM 
M & S simulation. 

The differences in geometry between the MLS and conventional ILS config- 
urations are associated with the location of the outer marker with respect to 
the runway touchdown zone (TDZ) . In theory, the final-approach common path 
can be shortened if RNAV-equipped aircraft utilize the wider, high-precision 
coverage afforded by the MLS system to intercept the ext ended- runway centerline 
closer to the threshold. Thus, for the TAATM M & S, the outer marker locations 
are 3.5 and 5.5 n.mi. from the TDZ for MLS and ILS, respectively, with corre- 
sponding differences in the locations of the approach GATE, ALTURA, NORTH2, and 
SOUTH2 waypoints. The altitude objectives at the gate (glide-slope intercept 
altitudes at the outer marker) are 1981 and 2134 m (6500 and 7000 ft) for glide 
slopes of 3.00^ and 2.75^ for MLS and ILS approaches, respectively. 

In general, the initial arrival clearance for the route segments from the 
perimeter to the lAF calls for a descent to the lAF altitude followed by a 
250-knot speed command. The delivery time for the speed command is computed 
to provide earliest time of arrival at the lAF based on perimeter penetration 
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However, if a holding stack is in use at the lAF, the next available stack alti- 
tude is assigned followed by a 21 0- knot speed conunand delivered at the appropri- 
ate computed time. All holding stacks consist of four levels starting at the 
lAF altitude and incremented by 305 m (1 000 ft) upward. 

To maintain the future controllability that would be lost by the 21 0-knot 
constraint in the holding pattern, a speed increase is usually requested by the 
speed-control routine (depending on existing conditions) upon departure from the 
lAF. This is an exception to the normal practice of issuing speed reductions 
only and is not objectionable from a pilotage viewpoint because of the potential- 
to-kinetic energy transfer available during the ensuing descent. 

The following sections describe the particular regions of application of 
the M & S control concept discussed in detail in section 3.0. Typical delay 
times attainable with speed control only are shown for a Boeing 737-class air- 
craft in the following list: 


BYERS arrival 66 sec 

LONGMONT arrival 61 sec 

SHAWNEE arrival 1 04 sec 

ELIZABETH arrival 65 sec 


2 . 1 BYERS Appr oach 


The straight-in BYERS approach navigation plate is shown in figure 2. Due 
to the limited speed-control delay capability and heavy utilization of this path 


FIXED PATH M&S/MLS 
DENVER, COLORADO 
BYERS Arrival - runway 26 

DENVER APR CON 120.5 288 1 

LOCAL CON 118.3 257.8 

ATIS 125.6 


N 

IVARIATION 13° EAST ~| 


DENVER 
1 16.3 DENI 




Cross at 
15QK 6,500’ 




7,000’ 


GATE 




ALTURA 


8 , 000 ' 



WATKINS 



Figure 2.- BYERS STAR navigation plate. 
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(40 percent of the arrival traffic}/ the aircraft on routes 1 and 2 are given 
landing priority over those on the other three approach paths. Thus, holding 
stacks at BYERS are required much less frequently than at the other lAP's. The 
remainder of the path contains two speed-control regions (BYERS to WATKINS and 
WATKINS to GATE) with firm sequencing and possible path stretching occurring at 
WATKINS. 


2.2 LONGMONT Approach 


Arrival routes 3 and 4, carrying 35 percent of the traffic/ merge at the 
MEEKER intersection before reaching the L(»JGMONT lAF as shown in figure 3. 


FIXED PATH M(ScS/MLS 

DENVER, COLORADO 

LONGMONT ARRIVAL - RUNWAY 26 

DENVER APR CON 120.5 288.1 

LOCAL CON 1 18.3 257.8 

ATIS 1 25.6 



Figure LONGMONT STAR navigation plate. 
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Speed cx>ntrol may be applied in two regions {LC»^GMONT to BRIGHTON and BRIGHTON 
to NORTH!), with possible path stretching at BRIGHTON. 

Firm sequencing takes place just prior to reaching NORTH! , at which point 
the minimum path direct to ALTURA is assigned if preceding traffic permits. 
Otherwise, the conventional downwind leg is assigned from which the initial DICE 
countdown to ALTURA is computed. Final DICE turns to GATE are computed from 
the minimum path or base (crosswind) leg along with a speed command such that 
final-approach speed reduction is completed at the gate. 


2.3 SHAWNEE Approach 

The lightest route loading (10 percent) occurs on arrival route 5 (SHAWNEE, 
fig. 4). Ironically, this path has the greatest speed-control capability, with 
three regions between SHAWNEE and SOUTH! . In addition, path stretching may be 
applied at CONIFER. Firm sequencing is done just prior to SOUTH! and employs 
the same options as described for the LONGMONT approach. 


FIXED PATH Mc^S/MLS 

DENVER, COLORADO 

SHAWNEE ARRIVAL - RUNWAY 26 

DENVER APR CON 120.8 363.0 

LOCAL CON 1 18.3 257.8 

ATIS 125.6 


N 


DENVER I 
1 16.3 DEN 



Figure 4.- SHAWNEE STAR navigation plate. 
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2.4 ELIZABETH Approach 


Arrival routes 6 and 1 , with 15 percent of the traffic, merge at the lAP 
(ELIZABETH, fig. 5) . Beyond ELIZABETH, the control regions are essentially the 
same as on the LONGMONT approach. 


FIXED PATH M&S/MLS 
DENVER^ COLORADO 
FI I7ABEYH arrival - 


RUNWAY 26 


DENVER APR CON 120.8 363.0 

LOCAL CON 1 18.3 257.8 

ATIS 1 25.6 


DENVER 
1 16.3 DEnI 


i 


Cross at 
150K 6.500’ 

^ ALTURA 


000 ' 


170K I 170K 

8 . 000 ' 8 . 000 ' 
SOUTH lT S0UTH2 


N 

IVARIATION 13° EAST I 



10,000- Afranktown 





KiOWA 
1 17.5 IOC 


Figure 5.- ELIZABETH STAR navigation plate- 


2,5 Controller Sector Allocation 

For the real-time M & S studies described in section 7-0, the Denver 
terminal area is divided into four controller sectors. The North Approach 
Control is responsible for aircraft on the BYERS and LCWGMONT paths between 
the perimeter and firm-sequencing points. Likewise, South Approach Control 
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is responsible for the same areas on the SPIAWNEE and ELIZABETH approaches. All 
control for the DICE areas from the firm-sequencing points of each route to 
the outer marker is allocated to Local Approach Control. Finally , the landing 
clearance and associated information is issued by Stapleton Tower Control. 

3.0 DESCRIPTION OF THE FIXED-PATH METERING 
AND SPACING (M & S) CONTROL CONCEPT 
3.1 Overall View of the Fixed-Path M & S System 

The basic functions of the M & S system are to establish a desired aircraft 
landing sequence, to assign scheduled times of arrival at the outer marker (OM) 
which provide adequate interarrival spacing between aircraft, and to offer con- 
trol commands in order to realize these scheduled times of arrival. In order 
to acconplish these functions the system utilizes three levels of control: 
Metering, delay spacing, and precise final spacing. The particular control 
algorithms used in this study are an adaptation of the procedures described in 
references 1 to 4 for the Denver terminal area. These algorithms are described 
in greater detail in sections 3.2 and 3.3. Flow charts that acconpany the dis- 
cussion in these sections can be found in appendix D. 

As seen in figure 6, the three levels of M & S control are performed on 
an arriving aircraft depending on its position relative to the runway. In the 
figure, the perimeter entry line represents the transition position from en 
route control to terminal control at terminal radar acquisition. The solid line 
from perimeter entry to the runway represents the desired fixed flight path of 
the aircraft within the terminal control region. The initial arrival fix (lAF) 
is denoted by a triangle on the flight path and is the point within the terminal 
area toward which en route control is responsible for feeding aircraft. The 
dashed lines in the figure represent the initiation of particular control actions 
which define the boundaries of the M & S control levels. Beginning at perimeter 
entry and moving toward the outer marker (OM) , the dotted lines represent the 
holding delay region, a backup delay vector region, and the fine-tuning direct- 
course-error (DICE) delay area, respectively. 

Metering is responsible for controlling the flow rate of aircraft into the 
terminal area. The metering region is bounded by the 20-min flight-plan-message 
and the 5-min flight-plan-message M & S control actions. This region falls 
within the en route control region where controllers are responsible for feeding 
aircraft to a terminal area lAF with adequate separation. The 20-min flight- 
plan message is issued by an en route controller to the M & S system when a par- 
ticular aircraft under the controller's jurisdiction is estimated to arrive at 
an lAF in 20 minutes. This message identifies the aircraft to the M & S system 
and triggers the metering control logic. This logic is responsible for calcu- 
lating needed en route delay when the flow rate toward the lAF or into the ter- 
minal area reaches some threshold above capacity. If no delay is requested by 
the system the en route controller permits the aircraft to proceed. If delay 
is needed, the M & S system issues to the en route controller a target time to 
deliver the aircraft at the lAF. The controller is then responsible for the 
aircraft achieving its target time. 
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Delay spacing is responsible for dynamically assigning and adjusting the 
landing sequence and scheduled outer-marker arrival times (SMT) and for issuing 
control commands in an attempt to match an aircraft's actual arrival time at 
the outer marker (OM) with its planned SMT. The delay-spacing region is bounded 
by the 5-min flight-plan message and the firm-sequencing M & S control actions. 
This region covers a portion of both the en route and terminal control regions. 
The 5-min flight-plan message from the en route controller to the M & S system 
represents an update of the estimated time of arrival (ETA) at the lAF for 
an aircraft. This update occurs after any required metering delay has been 
accomplished. 

Upon receiving a 5-min message for a particular aircraft, the M & S system 
assigns to the aircraft an initial scheduled outer-marker arrival time (SMT), 
which allows for adequate separation from previously scheduled traffic and which 
is based on a projected f irst-arrival-at-the-runway criteria. As the aircraft 
proceeds along its flight path, its ETA at the OM is periodically updated. 

Based on updated ETA's, the SMT's of all aircraft within the delay spacing 
region are dynamically adjusted in order to account for errors and maintain 
separation. This dynamic adjustment consists of scheduling modification and 
resequencing. 

At predefined points on the flight path within the delay-spacing region, 
attempts are made to match an aircraft's arrival time to its planned SMT by 
generating optional speed control, holding delay, and backup vector commands. 
Prior to the lAF, if a longer terminal delay is anticipated than can be accommo- 
dated by speed control alone, then a holding delay is used at the lAF. Between 
the lAF and the firm-sequencing point, delay is accomplished through speed con- 
trol with vectoring as a backup if longer delay is required. Speed control is 
accomplished by indicated-airspeed commands within a performance range for a 
particular aircraft type and are calculated for the estimated wind and the known 
remaining distance along the approach path. 

Nearer the runway, precise final spacing is used to correct for errors in 
aircraft performance during delay-spacing maneuvers in trying to meet SMT's. 

The firm-sequencing control action establishes a landing sequence which cannot 
be modified for aircraft within the precise final-spacing region. The SMT's 
within this region can only be modified by schedule slippage, which adjusts the 
spacing between aircraft by increasing their separation or closing unnecessary 
gaps. There are also normally two dynamic path-adjustment commands given in 
this region in order for an aircraft to achieve its SMT. These adjustments are 
accomplished using a direct-course-error (DICE) technique which determines the 
error at the OM which will occur if an aircraft immediately initiates a turn 
to a predetermined point on the path. 


3.2 En Route Functions 

3.2.1 General description .- En route control is responsible for feeding 
aircraft to the lAF's within the terminal area. The M & S system affects this 
control, however, by monitoring the flow rate of aircraft into the terminal area 
and limiting this flow, when necessary, by giving to the en route controllers 
target times for the aircraft to arrive at the lAF's. The en route controller 
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is responsible for maintaining separation between aircraft headed toward the 
same lAF/ identifying aircraft to the M & S system, and determining the delay 
maneuvers required to achieve any target times requested by the M & S system. 

The TAATM simulation of the M & S system includes a statistical modeling 
of the en route control region from the time an aircraft is identified to the 
M & S system until terminal area perimeter entry time. For each of the four 
IAF*s modeled for the Denver terminal area, an arrival queue of aircraft is 
maintained. Each queue is divided into two lists of aircraft, (1) aircraft for 
which 20-min flight-plan messages are pending and which have not yet been 
entered into the M & S control environment, and (2) aircraft which have been 
identified to the M & S system but for which the 5-min flight-plan message has 
not been delivered. Each list contains aircraft identification information, 
predicted control- function times at which an en route control action should be 
performed for each aircraft, and information concerning the control function 
desired for each aircraft. The control logic within TAATM is responsible for 
monitoring the queues, executing the desired function for an aircraft at the 
time specified within its queue, and placing aircraft in the appropriate queue 
list. Each queue represents a set of aircraft for which a single en route con- 
troller maintains responsibility. 

The traffic sample presented to the TAATM represents a desired set of air- 
craft to fly in the terminal area. The information associated with each air- 
craft in the sample includes its initial conditions (i.e., position, speed, and 
altitude), its performance type, desired route, and an offer time. For this 
study, the traffic sample was generated using a statistical model of arrival 
traffic in the Denver area. When the offer time for a given aircraft is reached, 
the model places the aircraft in the appropriate queue, dependent on the lAF 
toward which it will be flying. The control- function time associated with the 
aircraft in the queue is its offer time. The list within the queue, into which 
the aircraft is placed, contains aircraft waiting for a 20-min flight-plan- 
message control action. Aircraft for which a 20-min flight-plan-message control 
action has been successfully completed are moved within their queue into the 
list of aircraft waiting for a 5-min flight- plan-message control action. The 
control- function times associated with aircraft in this list are the predicted 
5-min flight-plan-message times calculated for the particular aircraft within 
the 20-min flight-plan control-action logic. Details of how these times are 
calculated are found in the following sections. 

The control logic which monitors the queues places highest priority on per- 
forming control actions for aircraft needing a 5-min flight-plan message. The 
assumption is made that the en route controller, responsible for all of the air- 
craft headed toward a particular lAF, would place highest priority on identify- 
ing to the M & S system those aircraft closest to the lAF. When the control- 
function time for a particular aircraft is equal to or less than the model's 
current clock time, the aircraft is flagged for the needed control action. If 
more than one aircraft within a queue needs a control action performed on a par- 
ticular simulated 4-sec radar scan, only the aircraft with the highest priority 
and also the earliest desired control-function time receives the control during 
that scan. Necessary control for other aircraft must therefore be delayed for 
at least one scan depending on the number of aircraft needing a control action 
performed during the given scan. 
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3.2.2 Twenty- minute flight-plan message .- The 20-min flight-plan message 
is defined as the message which identifies an aircraft to the M & S system when 
the aircraft is predicted to arrive at its lAF in 20 min. This message is 
delivered by an en route controller to the M & S system with the assumption that 
the controller has provided for adequate separation between this aircraft and 
the aircraft ahead flying to the same lAF. Since the offer times within the 
traffic sample are statistically generated and do not automatically insure this 
separationr the determination of the 20-min flight- plan-message time for a given 
aircraft requires simulated en route control logic to provide this separation. 

The aircraft having the earliest 20-min control-function time is considered 
by the control logic for a 20-min flight-plan message if there is no other air- 
craft in the same queue being considered for or receiving a 5-min flight-plan 
message. For an aircraft being considered for the 20-min message, the en route 
control logic determines if adequate separation at the particular lAF would 
exist between this aircraft and the previous aircraft accepted on the same route. 
If the aircraft cannot proceed and arrive at the fix conflict free, then the 
associated control- function time is rejected as a 20-min flight-plan-message 
time. Under this condition, the aircraft is left in its queue within the list 
of aircraft waiting for 20-min messages. The control-function times of all 
aircraft within this particular list, however, are then incremented by 4 sec 
(equivalent to one scan) and no further control is performed during this scan 
on any of the aircraft in this queue. If adequate separation exists between 
the aircraft pair, the control- function time associated with the aircraft is 
accepted as the 20-min flight-plan-message time T2 q* The M & S system then 
proceeds to determine the impact of the terminal traffic flow on this aircraft. 

The M & S system is responsible for calculating an estimated time of 
arrival (ETA) at the outer marker (OM) and for determining a metering scheduled 
marker time (MSMT) . For the purpose of metering the arrival of aircraft into 
the terminal area, the MSMT assigned to an aircraft is computed based on the 
MSMT assigned to the last aircraft accepted into the M & S system plus separa- 
tion time. This MSMT is used only to obtain a gross initial estimate of the 
flow of traffic and is not used by the system except during the 20-min flight- 
plan control action. No provision is made for a metering schedule or sequence 
update of the previous aircraft as it proceeds into the terminal area. 

The following equation is used to calculate the MSMT for the nth aircraft 
accepted into the system: 


MSMTn — max(MSMTj ^_1 + tg^p + + 15, ETAj^) 


The value is the required time separation at the OM between aircraft n 

and n-1 and is dependent on the performance and weight classes of the two air- 
craft. The value t^omp provides additional separation at the OM, if necessary, 
in order to compensate for the loss of separation due to speed differences along 
the common path on final approach. The 15-sec buffer is added to the required 
separation to compensate for the delivery errors of the system at the OM. For 
a normal error distribution and a separation violation rate of 5 percent, the 
value used would be 1.65 times the actual delivery interarr ival-error standard 
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deviation. However, the interar rival error is the very quantity which this 
study is measuring. Since the measure is not throughput rate, assuming some 
approximate value for the experiment yields a reasonable flow performance. The 
discussion in a later section on fast-time performance criteria indicates that 
interarrival-error standard deviation is independent of the time interval 
between scheduled arrival times. Since the buffer value only affects the sched- 
duled separation times, the interarrival-error standard deviation is independent 
of the buffer value assumed. The 1 5 sec used implies an assumed interarrival- 
error standard deviation of 9.1 sec, which the data show is reasonable. 

Metering the flow of aircraft into the terminal area is acccxnplished by 
calculating the required time delay t^^ in order for an aircraft to attain 
its MSMT and assigning an en route delay tg^ equal to if tj^ exceeds 

an input threshold E^: 


^eA “ ^rA 


(tRA > Ea) 


or 


^eA ” ^ 


(tRA = Ea) 


where 


tRA = MSMPn - ETAn 


If the value tgA greater than zero, the M & S system is responsible for 
requesting the appropriate en route controller to delay the aircraft by that 
amount. The M & S system does not, however, suggest a maneuver to achieve this 
delay. 

Completion of the 20-min flight-plan control action includes moving the 
aircraft within its queue to the list of aircraft waiting for the 5-min flight- 
plan message. The new control- function time associated with the aircraft is 
the time to deliver the 5-min message. The equation for cc«nputing this time 
is described in the next section. 

3.2.3 Five-minute flight- plan message .- The 5-min flight-plan message is 
defined as the message which identifies an aircraft to the M & S system when 
the aircraft is predicted to arrive at its lAF in 5 min. The time T 5 to 
deliver this message is computed within the 20 -min flight-plan logic and is 
stored as the control-function time in the appropriate queue. This time 
reflects any assigned en route delay, any errors which account for delivery 
accuracy between T 20 T 5 , and any necessary additional time to insure 

separation between aircraft pairs at the fix. The following equations are 
used to calculate T 5 : 
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T 5 = T 


(T > Tc ) 
^sep 


or 


^5 “ ^5sep ^sepA 


(T ^ 


Ts ) 

^sep 


where 


T = ^20 + + t20^ + tgA 

where ^ 20 ^ is a delivery error between T 20 and T 5 , "^Ssep predicted 

time at the lAF for the preceding aircraft on the same route plus separation 
time at the fix, ^sep^ i® ^ separation error at the fix, and 900 indicates 
a time value (15 min). The units in all of the preceding equations are seconds. 
The values of ^ 20 ^ and ^sep/^ random numbers based on a normal distribu- 

tion. The mean and standard deviation for ^ 20 ^ are, respectively, 0 and 
30 sec. The mean and standard deviation for ^sep^ are, respectively, 0 and 
2.5 sec . 

Simulation of the en route control region between the initial 20-min 
flight-plan message and the calculated 5-min flight-plan message time is 
restricted to maintaining the aircraft in its particular en route queue. The 
T 5 control-function times are maintained in en route queues depending on the 
lAF's toward which the aircraft are flying, and are ordered on increasing time 
within the appropriate list in the particular queue. No flight dynamics are 
simulated for an aircraft in the system until the terminal area perimeter entry 
time for the aircraft has been reached. This time is calculated during the 
5-min flight-plan message control action associated with T 5 . 

When the time for an aircraft to receive a 5-min flight-plan message has 
been reached, the model logic first determines whether or not there is room in 
the holding stack for the aircraft if it was to proceed to its lAF. If the 
stack is full, the aircraft is left in its en route queue list and all queue 
times within the list are incremented by 4 sec, the equivalent of a radar sweep 
(scan) time. This logic simulates the responsibility of an en route controller 
for maintaining a reasonable flow of aircraft to the lAF. If there is an avail- 
able altitude in the holding stack at the fix, the control-function time asso- 
ciated with the aircraft in its queue is accepted as T 5 and the control logic 
for a 5-min flight-plan message is executed. 

The 5-min flight-plan logic in the M & S simulation is responsible for 
updating the sequence of aircraft presently in the system and performing an 
initial scheduling of the new aircraft being accepted into the terminal area 
by the M & S system. The sequencing and scheduling algorithms used are the 
same as those described in section 3.3.2 for the terminal model. The only 
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delay option made available to the control logic at this timer however, is a 
hold at the initial fix. 

Transition of an aircraft from its en route queue to an active aircraft 
in the terminal area is accomplished by calculating a perimeter-entry time 
Tp, which corresponds to a radar-acquisition time for the terminal area. At 
time Tpr the aircraft is located at its initial terminal-entry point provided 
by the traffic sample, and flight dynamics are initiated which update the posi- 
tion of the aircraft at 4-sec intervals, simulating the radar scan time. The 
elapsed time required to fly from the perimeter to the fix Atpf is calculated 
based on flight-path data input to the model. Time Tp is then calculated by 
using the equat ion 


Tp = T5 + 


(300 - Atpf) 
300 


(300 + t 5 ^) 


where t 5 ^ is a delivery error between T 5 and the fix and 300 is a time 
value equal to 5 min. The value t 5 ^ is a normally distributed random number 
with a mean of 0 and a standard deviation of 1 5 sec. The elapsed time to 
fly between the perimeter and the fix is calculated by the ETA computation 
routine in TAATM described in reference 7. 


3.3 Terminal Functions 

3.3.1 General description .- The TAATM provides a realistic terminal area 
environment in which to test the M & S concept. In order to incorporate the 
M & S concept, modifications were made to the TAATM control logic and real- 
time simulated radar display. The RNAV and MLS navigational models were also 
added. The logic, which was retained basically intact from the original ver- 
sion of TAATM, includes the following: The flight dynamics of the aircraft 

whose positions are updated every 4 sec, the logic for assigning flight-path 
descriptions to an aircraft, the holding-stack maintenance logic, the communi- 
cation logic, the logic for determining when a control action should be per- 
formed for an aircraft, and the wind-model algorithms. Detailed descriptions 
of these logic components can be found in references 7 and 8 . 

The M & S concept incorporated in TAATM is a modification of the system 
in reference 4. Where possible, the logic described in reference 4 was adapted 
directly. Initial test runs, however, indicated that these procedures required 
modifications to handle problems generated by realistic traffic samples, holding 
procedures, flight dynamics, and errors which the TAATM model provides. The 
control philosophy implemented in TAATM follows the basic plan in reference 4 
of examining aircraft at discrete points within the terminal area in order to 
perform dynamic sequence and schedule updates and in order to determine whether 
a delay maneuver should be executed by an aircraft. Expansion of the sequencing 
concept, however, was necessary in order to avoid giving some aircraft unrealis- 
tic priority over holding aircraft. This modification will be described in 
detail later. The determination was also made that precalculated estimated times 


of arrival/ prestored delay times associated with given delay maneuvers, and 
fixed time separations at the OM cannot adequately handle arbitrary wind condi- 
tions and variability of arrival times because of the performance differences 
of the classes of commercial aircraft considered in our model. 

The M & S system establishes at fixed points, or at fixed times relative 
to fixed points, procedures which must be executed for a given aircraft as it 
proceeds to the runway. The basic philosophy assumes a single fixed point at 
the OM where the M & S system must maintain separation between aircraft. All 
control procedures prior to the OM must assist in maintaining this separation. 
The control procedures specified include updates of the landing sequence and 
scheduled outer-marker arrival times (SMT's) as well as dynamic adjustment of 
the speed and, if absolutely necessary, the path of the aircraft. Figure 7 
represents the LONCMONT approach to runway 26 at Denver and shows the types and 
relative positions of the control points which were implemented for this route. 
The perimeter-entry line represents the position at which the aircraft first 
appears to the tracking model following a ”hand-off’* from en route to terminal 
control at time Tp. The position of the aircraft is then updated every 4 sec 
(simulating a radar scan rate) as it moves along this path. Control actions 
are triggered as it proceeds to the runway. The terminal control elements in 
figure 7 are discussed in the following sections. 

3.3.2 Sequencing and scheduling .- The sequencing and scheduling algorithms 
of the M & S system incorporated in TAATM allow for the dynamic adjustment of 
the landing order and landing times of the aircraft in the terminal area as they 
proceed toward the OM. These adjustments are necessary when navigation, per- 
formance, surveillance, and wind errors cause the original assigned SMT*s to 
fall outside the capabilities of an aircraft on a given route. The schedule 
and sequence adjustment algorithms include the nominal schedule and sequence 
update algorithm, a resequencing algorithm, and a schedule slippage algorithm. 

An update of the landing sequence and SMT assignments occurs for all aircraft 
in the system which have not been firmly sequenced when any one of these air- 
craft passes a control action point that requires rescheduling. These control 
action points include the 5-min flight- plan message, radar acquisition (perim- 
eter entry), exit determinaticai for a hold, and firm sequencing. (See fig. 7, 
elements 2, 3, 5, and 7, respectively.) Resequencing of aircraft occurs when 
the normal sequence cannot be maintained because adequate separation does not 
exist between a pair of aircraft. Resequencing alters the normal sequence and 
permits an aircraft which has reached the limits of its delay capability to 
receive a higher landing priority. The system does not permit a change in the 
landing order of aircraft which have been firmly sequenced, but does permit 
these aircraft to initiate a fine tuning of their SMT*s in the form of forward 
or backward slippage. 

3. 3. 2.1 Scheduled outer-marker arrival time (SMT) computation: For each 

aircraft in the system from the 5-min flight-plan-message acquisition point to 
the runway, a set of parameters is stored in an array which describes the rela- 
tic«iship of that aircraft to the rest of the aircraft in the system, its delay 
capabilities in terms of bounds on possible times of arrival, and any delay 
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Figure 7.- Control actions performed on LONGMONT approach. 


which has already been assigned to the aircraft. The array is ordered on 
assigned SWTT's and at all times represents the latest update of the desired 
landing sequence and SMT assignments. The parameters stored in this SMT array 
which are used for SMT computation include the following: 


ETA, 


mi nr 


ETA 


max^ 


SMT^ 

Delay^ 


earliest estimated time of arrival at the OM of aircraft n, 
considering any delays assigned to the aircraft 

latest estimated time of arrival at the OM of aircraft n, 
considering the maximum delay capability which the system 
desires to use as an SMT bound 

assigned SMT for aircraft n 

required delay needed to realize the assigned SMT 

(Delayj^ = SMF^ - ETAj^, where ETA^ is the earliest estimated time 
of arrival at the OM excluding any hold delays assigned to the 
aircraft; if holds are assigned to aircraft n^ ETAj^ remains 
static until the aircraft is cleared out of the holding pattern) 


Initialization of an aircraft into the SMT array is acconplished during 
the 5-min flight-plan-message control action. The aircraft is simply placed 
at the end of the landing sequence and assigned an SMT which provides for ade- 
quate separation with the aircraft ahead in the sequence. (The separation cri- 
teria are discussed later in the section.) The SMT computation logic then 
redetermines the overall desired landing sequence and set of SNET's based on the 
delay capabilities of the aircraft within the array, including the newly entered 
aircraft. 


The procedure for determining the desired landing sequence is a modifica- 
tion of the logic described in reference 4. For the purpose of sequencing the 
aircraft, the terminal area is divided into three distinct regions instead of 
the two regions described in this reference. These three regions consist of 
the following: (1) All aircraft which have been accepted into the system from 

the 5-min flight-plan message to the holding fix, including all aircraft which 
have not yet received clearance past the fix; (2) aircraft which have received 
clearance past the fix but have not been firmly sequenced; and (3) all firmly 
sequenced aircraft. In reference 4, regions 1 and 2 are combined. 

On entry into the sequencing logic all aircraft are sorted into the three 
regions. Those aircraft which have not been cleared past the fix (region 1) 
are assigned a sequence ordered on the last ETA^^ which excludes assigned hold- 
ing delays. This establishes a f irst-ccme-f irst-ser ve sequencing order. For 
region 2 aircraft, the previously assigned sequence is maintained unless the 
control action being performed is the firm sequencing of an aircraft. If firm 
sequencing is anticipated for an aircraft, ETAj^^ij^^ is used to reorder up to 

three aircraft in region 2 behind that aircraft and all aircraft in front of 
it which have not been firmly sequenced. The order of those aircraft already 
firmly sequenced remains unchanged. 



The establishment of a desired landing sequence within regions 1 and 2 
prior to determining the overall landing sequence provides some stability 
between successive attempts to determine this sequence. If the separate order- 
ing within regions 1 and 2 is not first performedr two problems are generated: 

(1 ) The holding stacks create a discontinuity when attempting to order the air- 
craft on a projected first-come-first-serve (FCFS) principle. Aircraft on 
routes where no holding is being performed are then capable of receiving exces- 
sive priority over aircraft on other routes where holding is taking place, and 
this can create excessive holds for the holding aircraft. (2) The aircraft not 
yet to the fix are ordered on an FCFS principle based on ETAj^^j^^ which 

includes hold assignments. Under this condition, it is possible for aircraft 
headed for the same stack to have their landing sequence swapped back and forth, 
depending on when hold times are updated and the delay times assigned. 

The development of a desired overall landing sequence is accomplished by 
merging the aircraft from the ordered lists of aircraft from regions 1 and 2. 
While maintaining the order of these two lists, the ETAj^j^j^ value for each 
aircraft is used to merge the two lists on an FCFS principle. This enables the 
system to maintain a landing-sequence continuity within the regions while pro- 
viding a means to handle the merging of aircraft on a short route with those 
on a long route and the merging of faster aircraft with slower aircraft which 
are making a transition from region 1 to 2. 

Once the desired landing sequence has been developed, new SMT's are com- 
puted for all aircraft after the last firmly sequenced aircraft such that for 
an aircraft in order position n, 

SMTji ^ SMTn_i + SePn/n-1 


where 


ETAjnaXn = ^ ETA^iiin 

and Sepn/n«l is the time separation required between aircraft n-1 and n. 
The value Sepj^/i ^-1 is computed by 

^^Pn/n-1 ” ^sep ^comp ^ ^ 


where tg^p is the minimum required time separation at the OM between air- 
craft n and n-1 , is a separation at the OM to compensate for speed 

differences along the common path on final approach, and the additional 15-sec 
buffer is a separation to compensate for delivery errors of the system at the 
OM. These separation factors are identical to those used in en route metering 
as described in section 3.2.2. If an SMT cannot be assigned so that it does 
not violate the preceding criteria, then the desired landing order must be 
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itiodified in an attempt to develop a sequence and schedule which can be realized 
by all aircraft in the terminal area. This modification is referred to as 
resequencing. 

3.3. 2.2 Resequencing: Due to the limitations on the delay capability of 

some aircraft and dependent on the traffic interactions, the desired sequence 
of aircraft may not be realizable. If SMT's cannot be assigned under the cri- 
teria in section 3. 3. 2.1 for regions 1 and 2, resequencing is attempted. As 
SMT's are being assigned beginning with the aircraft behind the last firmly 
sequenced aircraft, the first aircraft encountered which cannot be assigned an 
SMT within its ETA bounds is swapped in sequence with the aircraft ahead of 
it. In other words, if aircraft n cannot fit behind aircraft n-1 , then 
the SMT for aircraft n-1 is ignored and an attempt is made to fit aircraft 

n behind aircraft n-2. If this second attempt is unsuccessful, then aircraft 
n continues to move ahead in sequence until it can be successfully assigned 
an SMT. Once a successful assignment is made, then the last aircraft which was 
passed in sequence by aircraft n becomes the next aircraft for which an SMT 
assignment is attempted. Priority is thus given to aircraft with limited delay 
capability. 

Under certain ETA boundary constraints, a resequencing attempt may still 
result in a failure to assign realizable SMT's to all aircraft under con- 
sideration. A failure mode is assumed during the resequencing process if more 
than six unsuccessful swaps are attempted or if no SMT can be successfully 
assigned to a given aircraft. Under either of these conditions, the updating 
of schedules is discontinued and the previously assigned SMT's are retained. 
Because at least one of these SMT's is unrealizable, the potential for a 
separation violation can then exist. Following each control action, however, 
a check is made to determine the magnitude in time of any potential inter arrival 
separation violation. This time error is displayed on the simulated radar 
display and the assumption is made that the problem will be resolved either by 
human intervention or by a later control action. 

3.3. 2. 3 Schedule slippage: For those aircraft which have been firmly 

sequenced, schedule slippage is initiated in order to fine tune the SMT's to 
account for minor errors as the aircraft are making their last turns to final 
approach and to close unnecessary gaps between aircraft. Schedule slippage is 
initiated prior to each direct-course-error (DICE) calculation performed for 

an aircraft. During schedule slippage no change in the landing sequence occurs 
and only those aircraft behind the aircraft initiating slippage are also con- 
sidered for slippage. 

Initially a desired SMT (DSMT) is calculated for the aircraft in sequence 
position n: 


DSMTn = max(SMTn_i + Sepn/n-1 / ETA^inj^) 

where SMTj^^-j is the SMT assigned to the aircraft ahead, Sepp/^-l time 
separation required between aircraft n and n-1 , and ETA^^^j^ is the earliest 
time of arrival for aircraft n. If 
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DSMTn = SMT^ 

where SMTj^ is the previously assigned SMT for aircraft n, then no slippage 
occurs. If 

DSMT^ > SMTj^ 


then backward slippage is required. Given that ETAj^^^ is the latest possible 
arrival time stored for aircraft n, if 

DSMTn ^ ETAjnaxn 

then the SMT can be and is slipped back to DSMT^^. If 
DSMT„ > 


then the SMT is slipped back to DSMTj^, even though it is not attainable. 

Beginning with the first aircraft under consideration for schedule slippage, 
each aircraft in turn behind this aircraft is tested to see if backward slippage 
is required. Backward slippage ceases when, for aircraft n, the DSMTy^ is 
less than or equal to the previously assigned SMTj^. If 


DSMTj^ < SMTj^ 


then forward slippage is possible. If 


DSMTn > ETAj„aXn 

then the aircraft cannot attain the DSMT^f but the SMT is still slipped forward 
to DSMTn- If 

DSMTn ^ ETA„,ax^ 

then a new SMT (SMTj^g^,) is assigned to aircraft n: 

SMTnew “ min(DSWrj^ + Slack, ETAj^g^j^, SMT^,) 
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The min function is used to prevent backward slippage if adding slack time 
to DSMT^ results in a greater time than either ETA^ax^ SMT^* The term 

Slack is a buffer (normally 10 sec) used during the direct--course-error (DICE) 
calculations on the backup-delay vector paths and on the downwind leg to prevent 
slipping an aircraft schedule forward and then having to slip it backward. This 
slack is removed during the final DICE turn on base. Each firmly sequenced air- 
craft behind the aircraft under consideration is tested in turn to see if forward 
slippage is possible for the aircraft. If an aircraft is encountered which is 
not in a DICE region or for which forward slippage is impossible, the testing 
of aircraft for forward slippage ceases. If forward slippage is possible for 
an aircraft behind the initially considered aircraft, a DICE calculation is 
triggered for that aircraft on the next 4-sec radar-scan update. 

3.3.3 Schedule maintenance .- Realization of a desired sequence and 
scheduled OM time is accomplished by determining necessary modifications to the 
flight profiles of the aircraft. Two types of schedule-maintenance logic are 
available within TAATM. Coarse schedule maintenance is performed by assigning 
speed changes, holding patterns, or backup-delay paths to the aircraft to 
achieve a desired SMT. Fine-tuning schedule maintenance is performed with cal- 
culated DICE maneuvers. 

The need for coarse schedule maintenance is determined at prespecified 
geometric points along the flight path or at calculated times during the flight 
for a particular aircraft. Control actions are performed which calculate esti- 
mated times of arrival at the OM and compare them to the aircraft's scheduled 
time of arrival. In addition, future delay capability for the aircraft is 
calculated and 70 percent of this capability is compared with the delay needed 
by the aircraft in order for it to achieve the scheduled time of arrival. No 
action is taken to modify the current flight-path segment if the delay needed 
can be realized at a future control-action point while retaining 30 percent of 
the future delay capability for subsequent accumulation of errors or schedule 
slippage. If the delay needed exceeds 70 percent of the future delay capa- 
bility, however, an action must be taken as soon as possible to modify the 
flight profile to achieve the necessary delay. 

Given that a delay action is necessary, an option is chosen from an array 
of possible options for coarse schedule maintenance associated with the par- 
ticular control action. Prior to the initial arrival fixes, no immediate delay 
option is available, but delay can be achieved by holding the aircraft at the 
fix. On the paths from the holding fixes to the final sequencing area the high- 
est priority option is a speed change. In this same region delay paths can be 
specified to achieve additional delay if a speed change alone is not adequate. 

The specification of options and the logic for determining which option to use, 
however, favors the utilizaton of speed-change delays over other delay manuevers. 
If there is no option available which is capable of achieving the needed delay, 
the aircraft is given the lowest speed profile and is permitted to proceed. 

If a conflict is produced because the delay problem cannot be resolved, a con- 
flict alert will be generated and displayed on the real-time simulated radar 
display as described in section 3.4. 
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The fine-tuning DICE calculation is handled as a separate control action 
type which is different from the coarse schedule-maintenance actions. Calcula- 
tions are performed to determine if an aircraft will arrive early or late at 
the QM if a turn to an RNAV waypoint is initiated immediately. Results of the 
DICE computation consist of a DICE value and a yes or no decision to initiate 
the turn. If the DICE value is 2 sec or less, an immediate/ direct-engage turn 
message is developed. This message is delivered by the communications logic 
to the aircraft during the next TAATM position update. DICE computations are 
performed in two basic regions. These are the path -stretching regions between 
the holding fix and the firm-sequencing points and the final-sequencing area 
between the firm-sequencing points and the gate for routes from the north or 
south. DICE is calculated in the path-stretching regions to determine when the 
aircraft should return to its nominal approach route if a delay path has been 
assigned to the aircraft. 

3. 3 .3.1 Speed option: When an arrival aircraft is entering a path segment 

with a speed-control option, the speed-control algorithm (SPDOPT) may be called 
upon by the control logic to determine the desired indicated airspeed for the 
aircraft to fly that segment. This determination involves calculating the esti- 
mated time of arrival at the outer marker and comparing it to the aircraft's 
scheduled outer-marker arrival time (SMT) . If the delay required for the air- 
craft to achieve its SMT is greater than 70 percent of the future delay capa- 
bility, then some speed modification will be applied on the current path seg- 
ment. Otherwise, action is deferred to a future control-action point. 

All the velocities relating to the aircraft and path-segment speed speci- 
fications are given in indicated airspeeds. However, in calculating times to 
fly, indicated-airspeed values are converted to equivalent ground speeds for 
the altitudes and wind speeds encountered by the aircraft. From this point in 
the speed-option discussion, the references to segment velocities will refer 
to the ground-speed equivalents. Once a segment ground speed is determined, 
then it is inversely transformed to desired indicated airspeed for the aircraft. 

Consider the events an aircraft goes through upon entering a single seg- 
ment of length S*| of an arrival route with initial ground speed . The 

arrival route is made up of an arbitrary number of n segments. The air- 
craft decelerates at a constant rate to some segment velocity which is 

within the speed range allowed for that aircraft type on the segment 
(Vjj^ = = Vl^ ) . For the balance of the segment in question the aircraft 

maintains the velocity Vg.^ . The purpose of the speed-control algorithm, if 

speed control is to be applied, is to determine the value of Vg^ which meets 

the necessary criteria. Once Vg^ is chosen, it then becomes the initial 
velocity for segment 2. The following sketch illustrates the segment-velocity 
relationship for a representative route with n segments. 

For the sake of discussion the initial assumption will be that there are 
no speed overlaps in successive segments, that is, Vl^ S more general 
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case will be covered later. Given , three estimated times of arrival at 

the outer marker are calculated external to SPDOPT for the value V3^ . These 
times are defined as follows: 

ETAR(V5^ ) outer-marker arrival time using V3^ on S-| and the highest 

velocities on the remaining segments 

XLLT(V3^ ) outer-marker arrival time using V3^ on S] and the lowest 

velocities on the remaining segments 

XMLLT outer-marker arrival time using the lowest velocities on all 

segments 

With the estimated times of arrival defined previously and V3^ = Vjj^ , the 
following can be determined: 

PDELAY present delay capability, XMLLT - ETAR(Vjj^ ) 

FDELAY future delay capability, XLLT(Vh^ ) - ETAR(Vh^ ) 

XLLT*(Vjj^) outer-marker arrival time by assuming no delay on S') 

and by using 70 percent of future delay capability, 

ETAR(Vh^) + 0 . 7 FDELAY 

The following diagram illustrates these quantities on a time-of-day axis. Also 
shown are the possible relative positions of these quantities to an assigned SMT 
for the aircraft. 
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There are four major cases which need to be considered depending on when 
the SMT occurs. These are as follows; 

Case 1 : SMT < ETAR(Vjj. ) . The aircraft is too slow to attain the assigned 

SMT/ even though the highest velocities are assumed. Rescheduling is therefore 
required to resolve the problem. (No speed-option delay is utilized for S-|.) 

Case 2: ETAR(Vh, ) ^ SMT ^ XLLT’ (Vh^ ) . In this case the aircraft will be 

assigned a Vpj^ to fly S-| since there is enough delay capability on the 
following segments to achieve the required delay. (No speed-option delay is 
utilized for S-|.) 

Case 3; XLLT’ (Vjj^ ) < SMT ^ XMLLT. In this case the speed needs to be 
determined in order to meet the SMT since there is not enough delay capability 
in the future. (Speed-option delay is necessary.) 

Case 4: XMLLT < SMT. The aircraft has insufficient delay capability to 

make its desired SMT. Under normal flight-error conditions/ some action should 
have occurred before reaching S*| to prevent this case. This case could occur/ 
however/ with the accumulation of severe flight errors. (Other delay options 
are examined to resolve the problem/ e.g./ backup-delay vectoring.) 

For case 3/ the criteria used to determine the segment velocty V 3 ^ used 
to meet the aircraft's SMT are: 


SMT = XLLT' (Vs^ ) 


and 


SMT = ETAR(Vg^) + k '[ FDELAY (Vg^ ) ] 
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where 


and 


FDEIAYCVs^ ) 


XLLT(Vg^ ) - ETAR(Vg^ ) 


k' = k + (1 - k) 


- Vs; 

-Vl, 


for k = 0.7. The speed-control algorithm determines Vg^ with the assumption 

that a portion of the delay required will be achieved on future segments. The 
value k defines the percentage (70 percent here) of the future delay capabil- 
ity which must be used in the future to help achieve the required delay if 
is assigned for S-| . The value k* is determined by a function which defines 
this percentage for the full range of possible speeds over S*| (i.e., Vjj^ to 

) . For the value of- k* is 0.7 and for Vl^ the value is 1.0. The 

inclusion of k* forces the highest reasonable Vg^ to be chosen and decreases 
the future delay capability held in reserve as Vg.^ approaches . 

In order to determine the value of Vg^ , it is necessary to solve the pre- 
ceding equations with XLLT(Vg^ ) and ETAR(Vg^ ) as explicit functions of 
Vg.j . These relations are contained in appendix B. The value Vg-^ obtained 
in this manner with the function k* will reserve some delay capability for 
schedule slippage. The previous time-of-day diagram only showed the relation- 
ship of ETAR, XLLT*, and XLLT for the value of Vg^ = . A more general 

representation as a function of Vg^ is shown in the following drawing: 



case 4 

case 3 

case 2 
case 1 
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Thus far in the discussion there has been no speed overlap on adjacent 
segments, that is, A more general procedure to handle the overlap 

situation (Vjj^ ^ = ^^2^ must be considered. It is possible that a 

Vs, can be determined such that Vs-| S The M & S system normally either 

maintains or reduces velocity for arrival aircraft. Whenever Vs^ ^ 
highest velocity permitted on segment 2 must be equal to Vs^ (Vjj^ = ^s-j / - 

Bounding affects ETAR(Vs^) since it is defined as the time to fly to 

the outer marker using Vs^ on S-| and the highest segment velocities on all 
remaining segments. Since the calculation for XLLT*(Vs^) involves using 

ETAR(Vs^), it is also affected. The relationships previously defined are ade- 
quate once the explicit function of XLLT(Vg^) and ETAR(Vg^ ) in terms of 

Vg^ for the overlap case are used. These relationships are contained in 
appendix B. 

3. 3. 3.2 Delay path stretching: The M & S system under evaluation in this 

report is based on a fixed-path speed-control concept with contingency delay 
path stretching. Prior to the aircraft reaching a holding fix, excessive delays 
which cannot be realized with speed control alone are handled by holding the 
aircraft at the fix. For aircraft which have passed the holding fixes, contin- 
gency delay path stretching is provided as an option to absorb excessive delay 
needed beyond the speed-control capability. 

Delay path stretching is accomplished by providing an alternate control- 
action description at specific points on the routes. In this case, the control 
logic returns to the main control routine and the new action is attempted frcxn 
the top of the control logic with a new path description and a new set of con- 
trol options. For the Denver M & S configuration this action is specified as 
an alternate path which provides a vector off the nominal path. A return to 
the nominal path is accomplished by DICE calculations after the aircraft has 
established itself on the vector. 

3. 3. 3. 3 Hold: Prior to an initial arrival fix (lAF) , the need for exces- 

sive delays can be resolved by holding aircraft at the fix. Holding stacks are 
provided at each of the four IAF*s so that up to four aircraft can execute hold- 
ing patterns simultaneously at each fix. The rules for holding an aircraft con- 
sist of the following: (1 ) Any aircraft which has previously been advised to 

anticipate a hold will in fact be required to hold, even if the traffic environ- 
ment changes prior to the aircraft reaching the fix; (2) an automatic hold is 
assigned to an aircraft if the aircraft ahead (assigned to same lAF) is holding 
or anticipating a hold; and (3) a hold is assigned if the delay required is in 
excess of 70 percent of the speed-control delay capability. 

Once the determination has been made that an aircraft will hold, a desired 
holding time is calculated based on the amount of delay needed minus 30 percent 
of the delay attainable with speed control in the future. This permits some 
forward slippage once the aircraft has left the hold. This desired holding time 
is then used to compute an actual assigned holding time. If the desired time 
is less than a minimum hold time, the minimum hold is assigned. 
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Hold assignments can be initiated by the 5-min flight-plan message^ at the 
radar-acquisition pointr or by the schedule-maintenance control action prior 
to the holding fix. Once the hold has been assigned, an update of the holding 
time can occur on subsequent control actions prior to the aircraft reaching the 
fix. Once the aircraft has initiated a hold, a holding-time update is possible 
only on the hold-exit-determination control action. This action occurs 50 sec 
prior to the aircraft's predicted time of departure from the holding fix. 

For the purpose of scheduling and sequencing aircraft, a possible holding 
delay time is included under certain conditions in the calculation of the maxi- 
mum estimated time of arrival. The inclusion of this delay time is dependent 
on the route being used by the aircraft and whether or not holds have already 
been assigned to the aircraft. For the straight-in BYERS route, this possible 
delay time is included only when the aircraft has already been assigned a hold. 
On all other routes, this time is included if the aircraft has not yet reached 
the holding fix, whether or not a hold has been assigned. This distinction is 
made between the BYERS route and all other routes in order to give BYERS air- 
craft a sequencing and scheduling priority. The result is that resequencing is 
attempted for a BYERS aircraft before a hold is considered. The reason for this 
priority is that the straight-in BYERS route has limited speed-control delay 
capability and is heavily used (40 percent of arrival traffic) . If a BYERS air- 
craft cannot be resequenced ahead of another aircraft, it may be forced to hold. 
Once it is assigned a hold, the route priority is removed by the inclusion of 
the maximum possible holding delay in the maximum ETA calculation of the air- 
craft. While BYERS aircraft are holding, all routes have equal priority. 

3. 3. 3.4 Direct-course-error (DICE): A DICE computation is performed to 

determine how late or how early an aircraft will arrive at the OM if a direct- 
engage turn to a new RNAV waypoint is initiated immediately. This direct- 
engage procedure uses the capability of a two-dimensional RNAV system (discussed 
in appendix C) to make a smooth turn from the present position and then fly 
directly to the next prespecified waypoint rather than continue toward the 
currently specified waypoint. A DICE control action can be initially triggered 
after an elapsed time following another control-action type or by a forward- 
slippage calculation performed on another aircraft. Once the calculation has 
been initiated the value returned from the DICE algorithm controls when the next 
computation is performed. The following conditions define how often the 
algorithm is called: 



Time before next 

DICE value, sec 

computation, sec 

00 

28 

^50 

8 

<50 

4 


The data used to determine these values was obtained from reference 4 and 
rounded to the nearest 4-sec radar scan. 
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The DICE computation can consist of an ETA calculation which assumes either 
a single immediate turn to a waypoint or an immediate turn towards a waypoint 
followed by a nominal second turn to the next waypoint. In either case the 
algorithm is responsible for determining the path to be flown, calculating the 
estimated time of arrival at the OM, and comparing this estimate with the sched- 
uled time of arrival (SMT) of the aircraft at the outer marker. Initial input 
variables to the algorithm consist of the following: 

Vt initial true airspeed 

Xq north radar coordinate for aircraft 

Yq east radar coordinate for aircraft 

hi initial track heading of aircraft 

Tj 4 time required to deliver direct engage message 

Tr roll-in time for turn 

All angles and headings used in this discussion are in degrees. In reference 
to figure 8, the predicted flight path of the aircraft is computed with the 
following equations: 


Di = Vi (Tm + Tr) 
e = mod (450 - , 360) 

X-i = Xq + Di sin 0 
Y-j = Yq + Di cos 0 


where is the distance traveled before the turn^ 0 is the angle measured 

counterclockwise frcm east to the heading vector, and (X*! ,Y-| ) is the initial 
position of the aircraft before the turn. If the turn direction Q] to way- 
point (Xp,Yp) is known, where Q] = T indicates a right-hand turn and Q] = -1 
indicates a left-hand turn, the position (X 2 /Y 2 ) at the end of the turn is cal- 
culated as follows: 


X2 = Xi + AXi 
Y2 = Yi + AYi 
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where 


Ax-j = RiQ-j (sin h 2 - sin h-| ) 

AYi = RiQi (cos h-| - cos h 2 ) 

The following equations relate to figure 8: 

R] = Vi2/(g tan $nqm) 

where 

g gravity 

nominal bank angle assumed for an RNAV direct-engage turn 


— X"! — R] sin h^ 

Yj^ = + R] Qi cos h'j 



01 = 90 - Qi (h-| - T-| ) - 

h2 = h] + 6i 0] 


If h 2 is negativef then h 2 = h 2 + 360. The distance traveled in the turn 
Dip^ is 

Dt^ = Ri^iCr 

where Cj^ is a conversion factor frcan degrees to radians. 
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If a single DICE turn is specified, the distance Dp to travel from the 
end of the turn (X2fY2) to the waypoint (Xp,Yp) is 

Dp = - Xp)2 + (Y 2 - Yp)2 

The total distance D-| to be flown from the aircraft's present position to the 
waypoint is 

D] = Dj^ + + Dp 

For a two-turn DICE computation, shown in figure 9 , the algorithm calcu- 
lates a point (X3,Y3) where the turn to intercept the nominal path to waypoint 
G begins. This calculation is accomplished by determining the offset to begin 
a turn from point (X2fY2> with initial heading h2 and terminating tangent to 
a line parallel to the nominal path to point G. The nominal-path heading to 
point G is denoted h3* By assuming a true airspeed V2f associated with 
the second turn, the radius of turn R2 is 

R 2 = V22/(g tan $nOm) 

If the turn direction Q2 to waypoint G is known, where Q2 = 1 indicates 
a right-hand turn and Q2 = -1 indicates a left-hand turn, the position (Xb,Y|3) 
is 


Xb = X2 + AX2 

and 


Yb = Y2 + ^^2 

where 


AX2 = R2Q2(sin h3 - sin h2> 

AY2 = R2Q2(c°® ^2 “ ^3) 
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Point (X^fYc) is the intersection of the line through (Xj^/Yj^) with heading 
h3 and the line through (X2fY2) with heading h3 + 90 . Point (XjjjY^j) is 


Xd = Xg + X2 - Xc 


Yd = Yg + ^2 - Yc 


The offset location (Xd^Y^) is used to define a line parallel to the nominal 
path to point G with heading h3 and intersecting the path to point (Xp,Yp) 
with heading h2* The intersection point is point (X3,Y3) . 

The distance traveled in the turn 0^2 is 


DT2 " R2®2Cr 


where 

02 = lAh'l 
02 = 360 - Ah' 
©2 = 360 + Ah' 

for 

Ah' = h3 - h2 
Point (X4,Y4) is 


(-1 80 ^ Ah' <1 80) 
(Ah' ^ 1 80) 
(Ah' < -180) 


X4 = X3 + AX2 
Y4 = Y3 + AY2 

The distance Dp flown on heading h2 is 
Dp = |/(X3 - X2)2 + (Y3 - Y2)2 
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The distance Dq flcwn from (X4,Y4) to point G is 


D(3 = |/(Xq - X4)2 + (Yq - Y4)2 

The total distance flown D“| from the aircraft's position to (X 3 ,Y 3 ) for the 
first turn is 


D] = + Dp 

The total distance flown D 2 from (X 3 fY 3 ) to point G for the second turn 
is 

^2 = ^T2 


These calculations provide the distances and headings necessary to calcu- 
late an estimated time of arrival (ETA) at the outer marker. A detailed 
explanation of ETA calculations in TAATM can be found in reference 7. The DICE 
value (in seconds) is calculated by DICE = SMT - ETA and indicates how early 
(if positive) or how late (if negative) the aircraft will arrive at the OM. 

The DICE value is then presented in the aircraft tag on the display. When this 
value is less than or equal to 2 sec, a message is generated and delivered to 
the aircraft to perform a direct engage to the next waypoint. 


3.4 System Output and Display 

COTimands and other information generated by the TAATM M & S control logic 
are output by alphanumeric display at the controller console in the form of tab- 
ular lists and aircraft tag components. It is envisioned that in a real system 
the controller would use the tag information to formulate the desired ATC com- 
mands. It should be noted that the display format adopted for TAATM does not 
necessarily replicate any existing or proposed actual air traffic controller's 
display but was developed for the intended uses of the TAATM program in an 
experimental envi ronment . 

The TAATM simulation has the feature of displaying abbreviated ccxtimand 
messages in ATC phraseology for a period of time dependent on message length. 
These messages and their delivery times are generated in accordance with the 
schedule-maintenance algorithms described in section 3.3.3. 

Two other types of lists are displayed as long as the information is 
applicable. One of these is a list for each holding stack of desired fix 
departure times of aircraft which have been assigned delays at the arrival 
(metering) fixes. The other is a conflict-alert list containing the aircraft 
flight numbers and predicted time errors for which conflicts are not currently 
resolved. This list and the data tags of the aircraft involved are displayed 
in a flashing mode for emphasis. 
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In close proximity to each aircraft position symbol, a data tag containing 
both conventional status and M & S control-directive information is written. 

The status data include the airline flight number, a heavy performance-class- 
aircraft indicator, a controller hand-off symbol, and radar altitude and ground 
speed. The control directives may include assigned holding times, a left- or 
right-turn vector instruction with a desired heading, a desired indicated- 
airspeed objective, and a DICE countdown value and an associated waypoint sym- 
bol and name. The name of the waypoint toward which the aircraft should be 
directed to turn as a DICE value approaches zero is displayed continuously when 
the value is less than 30 sec and flashing for values less than 1 5 sec. If the 
control logic determines that a waypoint should be entirely skipped by the air- 
craft, a leader line is flashed from the aircraft symbol to the desired way- 
point. 


Figure 10 shows a representative TAATM display and figure 11 describes the 
aircraft-tag nonenclature based on the data-block organization in reference 4. 


3.5 Interaction of the TCV B-737 Aft Flight Deck Simulator 
and the M & S Terminal-Area Simulation 

The TAAIM real-time simulation program and a fixed-base cockpit simulator 
form the basis of an overall mission simulator for real-time pi lot- in- the- loop 
experiments. The cockpit simulator used in the experiments is a replica of the 
aft flight deck of the TCV B-737-100 aircraft described in reference 9. The 
computer program for the aerodynamics, displays, and controls of the cockpit 
is programmed for the Control Data CYBER computer systems and Adage Graphics 
Systems. This simulation package includes a nonlinear mathematical model of 
the B-737-100 aircraft with the addition of landing-gear dynamics, gust/wind 
models, nonlinear actuator models, and instrument landing system (ILS) and 
microwave landing system (MLS) sensor models. 

Initialization of the cockpit takes place through a prescribed procedure 
in which the data for a selected TAATM aircraft is used to set the proper state 
values in the TCV program. Mode control for the TCV program is then slaved to 
the TAATM program for the ensuing operation in synchronized real time (SRT) . 

The TAATM program and the cockpit simulation are run simultaneously on separate 
Control Data CYBER 1 75 computer systems under the direct communication capa- 
bilities and constraints provided by the Langley Real-Time Simulation System. 
Reference 10 is a guide for the basic operation of the TAATM program and 
reference 11 contains documentation of the interface and operation procedure 
in the mission simulator mode. A block diagram of the real-time experimental 
configuration is shown in figure 12. Data transmitted between the programs 
consist of aircraft initialization parameters and mode-control signals from the 
TAATM and updated aircraft parameters from the cockpit simulation to TAATM 
during SRT operation. Data transfer of the updated aircraft parameters occurs 
at a rate of 32 times per second while operating. The TAATM, however, only 
examines the data lines at 4-sec intervals as if the inputs were from an 
automated radar terminal system (ARTS) radar. Tables 1 and 2 present the 
parameters involved in the mission-simulator data exchange. 
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TERMINAL AREA SIMULATION 


Gfir 18:46:46 



Figure 10.- EXcunple of TAATM simulated controller display. 
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Data-tag organization 




Data-tag 

fields 

Field 1 
Field 2 

DICE countdown 

Waypoint symbol 

associated with DICE 

(e.g., +75 ->■ 75 sec early at waypoint 
-5 ->■ 5 sec late at waypoint) 

(e.g., G ^ Waypoint at point G) 

Field 3 

Aircraft flight 
identification 


Field 4 

Symbol (H) for 
heavy aircraft 


Field 5 

A1 ti tude/100 


Field 6 

Hand-off symbol 


Field 7 

Ground speed/10 


Field 8 

Vector command 
Hold command 
Waypoint name 

(e.g., L210 ->■ Turn left to 210°) 
(e.g., HLD) 

(e.g., GATE ^ Turn to gate) 

Field 9 

Speed command 
Holding time 

(e.g., S210 ^ Speed 210 knots) 
(e.g., 04 -> Hold for 4 min) 


NOTE: Waypoint name in Field 8 can extend into Field 9 


Figure 11.- Aircraft data-tag information for TAATM display. 
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(CYBER 175) 
B-737 

simulation 


Instrument and control signals 
VOR/DME, MLS information 


DAC & ADC 


TCV B-737 simulation 


(CYDER 175) Position 
Terminal-area 


air traffic 
simulation 
model 


TAATM 

control 

console 


Control signals 


Position data 
and 

message information 


TAATM 

display 

console 


Audio link 


TAATM simulation 


Figure 12.- Block diagram of the real-time experimental configuration. 



The TAATM simulation treats the cjockpit as if it were part of the inter- 
nally generated traffic and will resolve conflicts and/or take advantage of 
situations created by controller /pilot actions. However, gross errors by either 
party can exceed the controllability of the system which will, of course, result 
in unresolved conflicts. 

In addition to the computer data transfer, voice communication is main- 
tained between the flight crew and the person acting in the role of the air 
traffic controller at the TAATM simulated radar display. During the mission 
simulation e^eriments, the TAATM controller and the cockpit crew interact in 
the manner described in section 3.6. 


TABLE 1 INITIAL ODNDITIONS TRANSMITTED FROM TAATM TO THE TCV SIMULATICW 


Data description 

Data range 

Resolution 

X-position, ft (m) 

±501 000 (±1 52 705) 

±0.06 (±0.02) 

Y-position, ft (m) 

±501 000 (±1 52 705) 

±0.06 (±0.02) 

Altitude, ft (m) 

0 to 30 000 (0 to 91 44) 

±0.90 (±0.27) 

Heading, deg 

0 to 360 

±0.01 08 

Indicated airspeed (IAS), knots . . . 

0 to 500 

±0.01 5 

RNAV bias errors: 



Azimuth, deg 

±5 

±3 X 1 0-4 

Range, n.mi 

1 

±1 

±6 X 1 0-5 


TABLE 2.- AIRPLANE -STATE PARAMETERS TRANSMITTED FROM TCV SIMULATION TO TAATM 


Data description 

Data range 

Resolution 

X-position, ft (m) 

Y-position, ft (m) 

Altitude, ft (m) 

Heading, deg 

IAS , knots 

Ground speed, knots 

Bank angle, deg 

±501 000 (±1 52 705) 
±501 000 (±1 52 705) 
0 to 30 000 (0 to 9144) 
0 to 360 
0 to 500 
0 to 500 
±90 

±0.06 (±0.02) 
±0.06 (±0.02) 
±0.90 (±0.27) 
±0.01 08 
±0.01 5 
±0.01 5 
±0.0054 
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3.6 Procedural Description of Human Interaction 

For the computer-assisted M & S system to operate smoothly, both the 
controller (s) and the pilot (s) must respond to the computer-generated requests 
within time ranges assumed by the ground computer. Although minor discrepancies 
are unavoidable due to human and mechanical factors, extreme navigation errors 
and/or ignored instructions can create circumstances to which the control logic 
may be unable to adjust. The model will then generate excessive resequencing 
and delays in its attempt to control the traffic. 

3.6.1 Controller responsibilities . - The primary responsibilities of the 

controller with respect to the computer-assisted M & S system are to interpret 
and relay the generated requests to the aircraft under his control, to monitor 
the traffic flow for safety, and to resort to alternative ATC procedures in the 
event of computer failure or catastrophic errors. In general, there are three 
types of control areas for each arrival route or STAR. These are: (1) Altitude, 

speed, and fix departure- time clearances (possibly including holding-stack 
delays) for the area between the terminal area perimeter and initial arrival 

fix (lAF) ; (2) fixed-path spjeed requests with possible delay vectoring between 
the lAF and firm-sequencing point; and (3) direct-engage and final-speed reduc- 
tion requests for fine tuning between the firm-sequencing point and the outer 
marker. The direct-engage requests to specified waypoints are issued in 
response to the computer-generated DICE-procedure countdown in the final control 
area and also for an imposed delay vector. Of course, other information such 
as traffic advisories, airport conditions, and radio-frequency changes (handoffs) 
is issued as required. 

3.6.2 Flight-crew responsibilities .- The nominal flight path, waypoint 
altitudes, and speeds are depicted on the navigation plate (STAR) for each 
approach path. An example procedure for the SHAWNEE STAR is shown in figure 13. 
The normal flight-crew procedure in flying an approach is to adhere to the STAR 
specifications except in response to deviations such as speed changes and holds 
or delay vectors requested by the controller (i.e., "control by exception"). 

To the best of their ability, the crew must determine the holding-pattern leg 
geometry necessary to meet the desired fix departure time for an assigned hold. 
Speed changes should be initiated immediately upon passing the points indicated 
on the approach plate or in response to a controller request. Altitude objec- 
tives should be met at or prior to passing the specified waypoint but generally 
have second priority to meeting speed objectives. 

Other time-critical responses include altitude changes within a holding 
stack, delay vector initiation, and direct-engage turns to waypoints in the 
DICE control areas. A direct-engage turn consists of simply breaking off from 
the current flight path and turning to a heading that will take the aircraft 
directly over the navigation waypoint specified in the controller request. 
Deceleration and descent rates should be consistent with those values modeled 
in the M & S computer for ETA computations. 
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DENVER, COLORADO 
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Figure 13.- Navigation plate for SHAWNEE STAR. 


4.0 ERROR MODELS 

The errors used in the study are a combination of the Terminal Area Air 
Traffic Model (TAATM) system and aircraft errors together with errors for the 
RNAV model supplied by the MITRE Corporation. The significant navigation models 
used for the M & S study are VOR/DME RNAV and MLS. The RNAV bias errors pre- 
sented in table 3 include angle and range errors associated with both the air- 
borne equipment (A6g and Ap^) and the ground station (A6g and APg) . They 
are combined in accordance with the equations 


^®FNAV - ^®g + ^®a 
^PRNAV - ^Pg + ^Pa 


where ^Srnav bearing-angle error and range error. 

The airborne-equipment errors are applied as a constant for a particular air- 


43 




craft and the ground- station errors are applied as a constant for a particular 
station site. The values and are applied to the exact posi- 

tion of the aircraft on a given scan in order to simulate the position measured 
by the airborne navigation system. 


TABLE 3.- RNAV BIAS-ERROR M3DEL 


Error 

Symbol 

^o value 

Application 

Airborne equipment: 




VOR 

A9a 

0.75° 

Determined once 
for each aircraft 
and kept constant 
during flight 

DME 

APa 

0.3 n.mi. 

Determined once 
for each aircraft 
and kept constant 
during flight 

Ground station: 


i 


VOR 

A6g 

1.0° 

Determined once 
for each ground 
station and kept 
constant 

DME 

Apg 

0.05 n.mi. 

Determined once 
for each ground 
station and kept 
constant 


The MLS errors presented in table 4 include signal noise (A0j^ and APj^) 
and bias errors (ABj^ and Apj^) which are applied depending on the range of the 
aircraft fron the station. The errors were chosen to model the basic wide- 
aperture/ high-resolution MLS system intended for high-density terminal areas. 
For this system two MLS coverage configurations, ±40^ and ±60^ coverage, were 
used in the experiments. Within this coverage the range errors for both bias 
and noise vary as a function of distance frc«n the MLS station site. The errors 
are combined in accordance with the equations 

^®MLS “ ^®b 

APmls = Apb + APn 

where A6j^S is the angle error and Apj^^S is the range error. The values 
A0MLS Ap|yiLs are applied to the exact position of the aircraft on a 

given scan in order to simulate the position measured by the airborne navigation 
system. 
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TABLE 4.- MLS ERROR MDDEL 


Error 

Symbol 

la value 

Application 

Bias: 




Azimuth 

A6b 

0.066° 

Determined once for 
each aircraft and 
kept constant 
during flight 

Range 

APb 

4.3 m 
(46.8 ft) 

t>88.7 m 
(291 .0 ft) 

Determined once for 
each range for 
each aircraft and 
kept constant for 
the balance of that 
range 

Noise: 




Azimuth 

c 

CD 

< 

0. 091 o 

Redetermined every 
4 sec for each 
aircraft 

Range 

APn 

^3.7 m 

(12.2 ft) 

bg.8 m 
(22.3 ft) 

1 Redetermined every 

4 sec for each air- 
craft and applied 
dependent on range of 
aircraft from station 


^For range ^4.5 n.mi. 
^or range >4.5 n.mi. 


The aircraf t-per formance deviation errors presented in table 5 are used 
to simulate the precise achieved values as opposed to the exact desired values. 
The first four deviation errors in the table are redetermined for each aircraft 
whenever a speed or altitude change is initiated. The application of speed and 
altitude deviation is discontinued once the aircraft is on final approach. The 
application of a descent-rate deviation is discontinued once the aircraft is 
between the outer marker and the runway. The last deviation included in the 
table is the time-delay flight technical error. This error simulates pilot 
variability in initiating a DICE turn. The assumption is that the pilot will 
respond to a direct-engage message from the controller with a mean time of 4 sec 
and a standard deviation of 3.5 sec. 
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TABLE 5.- AIRCRAFT-PERFORMANCE DEVIATION ERRORS 


Error 

la value 

Application 

Altitude, m (ft) 

30.5 
(1 00.0) 

Determined each time 
altitude change is 
initiated 

Speed, knots 

5 

Determined each time speed 
change is initiated 

Descent rate, percent 

5 

Determined each time a 
descent is initiated 

Deceleration rate, knots/min . . . 

5 

Determined each time a 

deceleration is initiated 

Flight technical time deviation 
in initiating a turn, sec . . . 

3.5 

Determined each time a 
DICE turn is calculated 


The wind-speed and direction bias errors are applied as straight percent- 
ages of the wind values which affect an aircraft's position on a given scan* 
The wind model and its associated errors consist of the following: 


Wg = [w + AsgCh - ho)] (1 + Bg) 

Wd = ['1' + Y(h - ho)] (1 + B<3) 


where 

Wg actual wind speed 

0) assumed wind speed at ground level 

assumed change in speed with change in altitude 
h aircraft altitude, MSL 

hQ airport altitude, MSL 

Bg wind-speed bias 

Wj) actual wind direction, deg true 

^ assumed wind direction at ground level, deg true 
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Y 


assumed change in wind direction with change in altitude 


wind-direction bias 

The ncmiinal values for Denver consist of the following: 

03 = 7 .91 8 knots 

ASg = 0.007769 knots/m (0.002368 knots/ft) 

y = 2770 

Y =0 deg/m (0 deg/ft) 

ho = 1 625 m (5330 ft) 

Bg =10 percent 

B^ =5 percent 

The en route metering errors presented in table 6 are used in the determi- 
nation of the 20-min and 5-min flight- plan-message times. The first two errors 
correspond to the uncertainties in the estimated times of arrival. When a min- 
imum separation is called for between the 20-min and 5-min flight-plan messages, 
the en route separation- time error corresponds to a variation in the inter- 
arrival separation time conpared to the minimum required. 

TABLE 6.- EN ROUTE METERING ERRORS 


Error^ 

la value, sec 

Arrival-time error 

30 

20 min from fix 


Arrival-time error 

1 5 

5 min from fix 


En route entrail 

in 

CM 

separation error 


^Determined once for 

each aircraft. 


The surveillance errors presented in table 7 are used to simulate the 
difference between the actual position of an aircraft and the radar measurement 
of the position. The error values are redetermined for each aircraft on each 
simulated radar scan . 
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The determination of an aircraft's actual final- approach and touchdown 
speeds is based on the aircraft's gross landing weight. Since weight is related 
to the load factor of the aircraft, there exists a functional relation between 
final speeds and load factor. By using this relation, the TAATM aircraft final- 
approach and touchdown speeds are determined using a mean load factor of 50 per- 
cent and a standard deviation of 1 5 percent. The load factor for a particular 
aircraft is determined based on a normal distribution which is bounded by 1 0 
and 1 00 percent. If a load factor is computed outside these bounds a new load 
factor is recomputed. The factor chosen for a particular aircraft remains con- 
stant over the duration of its flight. For the purpose of estimating times of 
arrival, the control logic assumes a constant 50-percent load factor for all 
aircraft. 


TABLE 7.- SURVEILLANCE ERRORS 


Error^ 

1 a value 

Range 

79.2 m (260.0 ft) 

Azimuth 

0.25° 

Mode C altitude 

1 2.5 m (41 .0 ft) 

i - 


^Determined each scan for each aircraft 
position update. 


5.0 IDENTIFIED LIMITATIONS AND CONSTRAINTS OF THE 
TESTED FIXED-PATH M & S SYSTEM 

During the M & S study several limitations and constraints on the system 
were identified which are functions of the terminal geometry used, the proce- 
dures for separating aircraft, and the sequencing and scheduling control algo- 
rithm. The following discussion analyzes the effects of these limitations and 
constraints on the operation of the experimentally tested M & S system. Within 
the bounds of this study, however, no attempt was made to remove them from the 
system. 


5.1 Geometry Constraints 

The terminal geometry described in section 2.0 includes the straight-in 
BYERS approach route. Unlike the remaining routes, the controllability 
available on this route does not include the two-turn DICE prior to final 
approach. The two- turn DICE provides the advantages of path shortening and 
path stretching to adjust to an assigned SMT close-in to the OM. It is also 
capable of accommodating SMT schedule slippages. In contrast, the control- 
lability on the BYERS approach includes a path-stretching vector which is only 
used as a back-up to its speed-control capability. The decision to initiate 
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the vector occurs only once on the approach route, at the final speed-control 
point. If an aircraft passes this speed-control point without initiating a 
vector, further controllability is not available for the aircraft and schedule 
slippages cannot be acconunodated. To lessen the impact of this constraint, 
the BYERS approach is handled as a priority route. For this route, hold assign- 
ments are applied prior to the fix only after resequencing has failed to resolve 
a separation problem. For all other routes, hold assignments have priority over 
resequencing. Additionally, the en route delay threshold for this route is set 
at 60 sec instead of the nominal 240 sec applied to the other routes. The 
impact of these differences in handling BYERS' aircraft is to smooth the flow 
of aircraft on this approach by forcing larger and more frequent en route delays 
and by minimizing holds. 


5.2 Separation-Procedure Limits 

The limitations and constraints with respect to the procedures for separat- 
ing aircraft consist of the following; (1 ) Separation is insured only on ini- 
tial arrival at the lAF's and on final approach; (2) the schedule- slippage pro- 
cedure can shift a separation problem from an aircraft which has generated the 
need for backward slippage to an aircraft behind with limited delay capability; 
(3) the ground rules for the assignment of holds can generate a cascade of holds 
resulting in excessive separations; and (4) the accumulation of aircraft flight 
errors can generate a schedule-maintenance problem which cannot be resolved by 
the system. The result is that not all separation violations and gaps can be 
resolved by the M & S system studied. 

5.2.1 Common-path separation violations .- Because separation is insured 
only on an aircraft's initial arrival at the lAF and then again on final 
approach, separation violations due to speed differences can occur between a 
pair of aircraft on a common path. These violations are possible but are rare 
because they occur only under certain circumstances for aircraft pairs exiting 
from holding patterns and for aircraft pairs with large speed errors. Because 
separation is insured on initial arrival at an lAF and on final approach and 
because the speeds of the aircraft are decreasing as they approach the OM, 
aircraft normally have excess separation upon leaving the lAF. The separation 
distance then closes gradually to the minimum required separation distance at 
the OM due to decreases in speeds. 

If a pair of aircraft are holding at the fix, it is possible for the 
schedule-maintenance logic to release the first aircraft assuming no speed- 
control delay will be used in the future to realize the SMT. The second air- 
craft can then be released from the hold assuming speed-control delay will be 
required. The result is a loss of separation just past the fix until the second 
aircraft receives its speed-control delay so that the separation distance is 
increased. 

Consider, for example, the LCWGMONT approach. (See fig. 3.) Assume that 
a minimum required 3-n.mi. separation distance between a pair of aircraft at 
the OM results in a time separation of 82 sec. Assume also that the delay 
achievable by using speed control is 62 sec (i.e., the speed delay for a B-737 
on LONGMONT approach) . The schedule-maintenance logic for hold-exit determi- 
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nation will permit the second aircraft to proceed past the fix if its needed 
delay can be achieved by using 70 percent or less of its speed-control delay 
capability in the future. In other words, the second aircraft can exit the fix 
needing to delay for 43 sec in order to meet its SMT. If the first aircraft 
exited fran the hold needing no delay, then the pair of aircraft is separated 
at the fix by approximately 39 sec. If both aircraft leave the fix at an 
indicated airspeed of 210 knots and an altitude of 4267 m (1 4 000 ft), the 
ground speed of each aircraft is approximately 257 knots, assuming no winds 
are present. With a separation time of 39 sec, the separation distance would 
be about 2.8 n.mi. The condition can therefore exist where the minimum 
3.0-n.mi. separation is not maintained over the entire route. 

A loss of separation between a pair of aircraft can also occur on a route 
because of large speed errors. The speed-error model computes actual speeds 
attained by aircraft which can be faster or slower than the M & S requested 
speeds. The possibility exists that a pair of aircraft on the same approach 
route have large speed errors and are scheduled to land consecutively. Depend- 
ing on the magnitude of the errors, and if the speed of the first aircraft is 
slower than desired and the speed of the second aircraft is faster than desired, 
a loss of separation can occur before the control logic is capable of recogniz- 
ing and resolving the problem. 

5.2.2 Schedule- slippage-procedure limitations . - Another type of problem 
can exist due to the schedule-slippage procedure. This procedure is designed 
to close gaps between aircraft and resolve possible separation violations at 
the OM by adjusting SMT*s without performing any resequencing of aircraft. 

This adjustment can consist of either a forward or backward slippage of the 
SMT's to resolve accumulated flight errors. Assume that an aircraft has accu- 
mulated flight errors such that SMT backward slippage is required because of 
the aircraft’s late estimated time of arrival at the OM. Once backward slip- 
page has been accomplished for this aircraft, all aircraft behind must be 
examined in turn and schedules must be slipped backward until slippage is no 
longer necessary or cannot be realized for a particular aircraft. If the 
second case occurs, the possibility can exist that an aircraft behind is penal- 
ized because of errors accumulated by the aircraft initiating the slippage. 
Separation violations occur then between a pair of aircraft which did not gen- 
erate the flight error requiring the slippage. This occurs particularly when 
one of the aircraft behind is on the BYERS approach and has already used up its 
limited delay capability. 

5.2.3 Holding-procedure constraints .- In contrast to the separation prob- 
lems mentioned previously, the holding logic is capable of generating separation 
gaps in the traffic flow on a single approach route. This occurs because of 
the severe constraints of the ground rules associated with the assignment of 
holds. As stated previously, the ground rules consist of the following: 

1 . Any aircraft which has previously been advised to anticipate a hold will 
in fact be required to hold, even if the traffic environment changes prior to 
the aircraft reaching the fix. 
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2. Hold assignments can occur during the 5-min flight-plan message action, 
at radar acquisition, at the control action just prior to the holding fix, and 
during the hold-exit determination. 

3. An automatic hold is assigned to an aircraft if the aircraft ahead 
(assigned to the same holding fix) is holding or anticipating a hold. 

4. A hold is assigned if the delay required is in excess of 70 percent of 
the speed-control delay capability. 

Values for 70 percent of the speed-control delay capability for the four 
arrival routes for a B-737 are presented in the following list: 


BYERS arrival "^46 sec 

LONGMONT arrival '-43 sec 

SHAWNEE arrival ^73 sec 

ELIZABETH arrival '-46 sec 


Therefore, if a B-737 on the LONGMONT approach requires a delay of greater than 
43 sec prior to the holding fix, the aircraft will be required to hold. A 
single hold for one aircraft can then generate a cascade of holds for all air- 
craft behind it which are proceeding to the fix. Resolution of the holds 
requires a sufficient break in the traffic flow to the fix such that the holding 
stack is cleared before another aircraft is offered into the active traffic for 
that same fix. A minimum hold for a B-737 requires approximately 3 min depend- 
ing on the altitude at which the hold is executed. If the aircraft requires 
a delay of 50 sec, an excess delay of approximately 1 30 sec is actually achieved. 
It should be noted, however, that resequencing can occur such that the gaps can 
be used by aircraft on other approach routes. 

5.2.4 Limited-delay-option specification . - Another limitation exists within 
the experimentally implemented M & S system because the specified delay options 
are finite. Depending on the magnitude of the flight errors, the traffic load, 
and the delay capability, the possibility exists that an aircraft at a control- 
action point after the holding fix does not have adequate delay capability 
remaining to achieve its SMT. The resulting individual-aircraft error may or 
may not, however, result in separation error depending on the relative error 
between aircraft pairs. Following each control action performed, a separation 
error predicted at the OM is calculated for all aircraft pairs within the termi- 
nal area. If an aircraft's SMT is either early with respect to its calculated 
minimum ETA or late with respect to its calculated maximum ETA, a potential 
separation error for the aircraft exists. By comparing the potential separation 
errors between an aircraft pair, a predicted separation error can be calculated. 
If the predicted error indicates a potential loss of separation, the desired 
aircraft delay to resolve the problem is displayed on the simulated radar dis- 
play. The simulation run, however, proceeds and the assumption is made that 
the problem will be resolved by human intervention or by M & S resequencing or 
rescheduling in the future. 
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5.3 Limitations of the Scheduling Algorithms 

Two important limitations have been identified within the control algo- 
rithms concerning scheduling. The first is that the knowledge of the delay 
capability used by the scheduling and sequencing logic is limited to a single 
pair of minimum and maximum anticipated ETA*s at the OM for a given aircraft. 

The second is that the en route metering function is performed independent of 
the terminal sequencing and scheduling algorithm. 

The first control limitation is a result of the types of ETA information 
made available to the scheduling algorithm. The algorithm for assigning an SMT 
attempts to bound the SMT by a minimum and maxium possible ETA at the OM for 
a given aircraft. These bounds represent the SMT limitations which are assumed 
at a given control-action point for the realization of an SMT. The time dif- 
ference between the two ETA's, however, does not necessarily reflect the total 
controllability available for the aircraft on the entire path. No additional 
information is presented to the scheduling algorithm that describes achievable 
ETA'S with path shortening or stretching. The result is that sequencing and 
resequencing decisions are made for the aircraft without a complete knowledge 
of the controllability of the aircraft. If SMT assignments cannot be made 
within the nominal ETA bounds, this additional information could be used to more 
realistically resolve the problem. 

Within the constraints of a fixed-path speed-control M & S philosophy, this 
scheduling procedure places highest priority on the realization of SMT's by 
speed control with a secondary priority on holding aircraft and resequencing 
SMT's. Vectoring is used by the schedule-maintenance algorithm only as a backup 
option to realize an SMT, but the scheduler does not assign SMT's with the 
assumption that vectoring will be used. If the scheduler was given more know- 
ledge of the controllability which could be achieved by path shortening and 
stretching, speed control could still be given highest priority but the second- 
ary priority would be placed on path changes rather than holding and resequenc- 
ing. With this change, the problems generated by holding aircraft could be 
reduced, but more vectoring would occur within the terminal area. 

A severe constraint of the studied M & S control philosophy is that the 
scheduling algorithm for performing en route metering is not coupled to the 
terminal scheduling algorithm. The en route metering function is designed to 
regulate the flow of aircraft into the terminal area but does not recognize the 
dynamically changing environment within the terminal area. En route delays 
are assigned based on a scheduling algorithm which differs from the terminal 
scheduler in the following two ways: (1) The terminal scheduler performs 

sequencing based on a projected f irst-come-f irst-serve-at-the-runway principle 
whereas the en route scheduler simply assigns a first-in-first-out order, and 
(2) the en route scheduler does not use the dynamically updated ETA information 
available to the terminal scheduler. The result is that large delays in the 
form of holds or vectoring are still necessary within the terminal area even 
though en route delays have been assigned. If, however, the en route scheduling 
algorithm was coupled with the terminal algorithm, potential problems within 
the terminal area could be anticipated more accurately in the en route area and 
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en route delays could be assigned accordingly. Gross delays could therefore be 
absorbed at the higher, more fuel-efficient altitudes. Results discussed in 
section 6.2 support this view. 


6.0 FAST-TIME STUDIES 

The primary goals of this study were to evaluate a particular fixed-path 
metering and spacing (M & S) system, described in section 2.0, and to determine 
if the M & S*s delivery performance was improved by using the microwave landing 
system (MLS) instead of the VOR/DME as the RNAV data source during approach, 
fine-tuning, and DICE directed maneuvers. Both ±40^ and ±60^ MLS were evaluated. 
Each aircraft was assumed to be equipped with a skip-waypoint, two-dimensional 
RNAV system operated as discussed in sections 3. 3. 3. 4 and 3.6.2. Data collected 
from fast-time (batch) runs of the computerized Terminal Area Air Traffic Model 
(TAATM) were used as a basis for evaluation. The data are analyzed and plotted 
in this section. 


6.1 Performance Criterion and Data Combination 

The aircraft inter arrival-error standard deviation at the outer marker is 
the primary criterion of performance evaluation and indicates the variation in 
interarrival error achieved by a system. Interarrival error is the difference 
between realized aircraft interarrival time and that of the expected interar- 
rival time obtained from the scheduled times between succeeding aircraft pairs. 
(See eq- (All) in appendix A.) This measure translates the errors of the M & S 
system in meeting its time objective into an indication of system separation 
performance, which is more closely related to ATC safety than the actual arrival 
errors. For instance, if all aircraft were to arrive with sizable time errors, 
but of the same magnitude and direction, then interarrival separation is still 
maintained. Another attractive feature of the measure is that it is theoret- 
ically independent of arrival rate. Time errors are determined from each air- 
craft's own specific schedule time to which the system is working to meet. 

Within limits, as the arrival rate increases, the interval between scheduled 
times decreases but the dynamics of meeting each aircraft's time objective 
remains the same. Therefore, results of separate data runs of the same system, 
even with different arrival rates, can be combined to obtain an improved pooled 
estimate of the variance of the system under study. 

Arrival-only traffic was used to increase the number of interarrival times 
for data analysis purposes. An estimate of the system inter arrival-error stand- 
ard deviation was obtained from a set of data runs with all parameters except 
arrival rate held constant. A data set consisted of 2 separate runs for each 
of the arrival-rate samples of 25, 30, 35, and 40 aircraft per hour. Each of 
the 8 runs in a data set contained 2 hr of steady-state, outer-marker inter- 
arrival data. A combined pooled estimate for the data-set variance was used 
to measure the system inter arrival-error variance. If k samples from a popu- 
lation having a common variance are available for estimating the variance, the 
pooled estimate Sq^ is defined by the following: 
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where is the sample-run variance and nj^ is the number of data points 

for the kth sample. The pooled estimate of standard deviation is simply Sq. 


6.2 Approach Navigation 

6.2.1 G eneral M & S performance .- Figure 14 shows the fast-time values of 
the interarrival-error standard deviation for various final-approach navigation 
systems. The figure shows a standard deviation in the range of 7.5 to 1 0 sec 
depending on the particular approach navigation used. Current manual ATC pro- 
cedures normally achieve values in the order of 1 8 to 21 sec (ref. 12). Thus 
the relatively fixed-path, speed-control M & S system studied gave reduced 
instrument flight rule (IFR) delivery dispersions in comparison with today's 
manual vectoring and yielded results comparable to M & S systems employing only 
vectoring. 

6.2.2 Total-geometry performance .- For the terminal route geometry studied, 
figure 14 shows comparison of the overall delivery performance of the M & S 
system using various approach navigation systems. There is no statistically 
significant difference between the performance of the ±40^ MLS system and that 
of the ±60^ MLS system. Use of the equal-tails F-test at the 1 -percent level 

of significance indicates there is a difference in the delivery performances 
between using VOR/DME and ILS and using MLS as the final-approach navigation 
system. Although statistically different, the standard-deviation decrease 
from 9.2 to 7.9 sec shows that MLS offers only a slight improvement in system 
delivery precision over VOR/DME and ILS for the particular M & S and terminal 
route geometry used. 

6.2.3 Route effect .- The total-geometry system performance discussed pre- 
viously is due to the interaction of two types of approach routes. It is 
desirable to separate and assess the individual performances of the straight-in 
routes as well as the other three routes, all of which use two DICE turns to 
supplement the time-control-by-speed instructions. Figures 15 and 16 illus- 
trate two types of traffic samples which were run with that end in mind. One 
sample contains only straight- in traffic and the other contains no straight- in 
traffic. Figure 1 7 indicates the inter arrival-error standard deviation delivery 
performance for the final- approach navigation system evaluated for the two types 
of traffic samples. The F-test at the 1 -percent level of significance shows 
that the two-DICE-turn route is worse than that for the all straight- in traffic 
when the VOR/DME and ILS system is used for approach navigation. When the 

MLS system is used, the DICE-turn routes yield a reduction of about 1.5 sec in 
the interarrival-error standard deviation compared with that from the all 
straight- in route. 

Figure 1 7 also indicates that if there had not been a straight-in route 
in the geometry but instead only routes using two DICE turns, then the VOR/DME 
and ILS and MLS delivery-performance difference would be larger than that for 
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Figure 14.- Interarrival delivery precision for various final- approach 
navigation systems. Wind error of 10-percent strength and 5-percent 
direction. Confidence interval of 95 percent for values shown. 
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Figure 15.- Traffic flow for 30 aircraft per hour data sample with all 
traffic on straight-in route. Final-approach navigation with 
±60° MLS. 
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Figure 16.- Traffic flow for 30 aircraft per hour data sample with 
all traffic on DICE-turn routes. Final approach navigation 
with ±60° MLS. 
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Figure 17.- Route delivery characteristics. Interarrival delivery precision 
for all-DICE-turn traffic and for all-straight-in traffic. No wind error; 
95 percent confidence interval for values shown. 
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the total-geometry case* With all-DICE-turn routes, the VOR/DME and ILS 
inter arrival-error standard deviation is 8.5 sec, whereas the MLS is 5.9 sec. 

So for the all-DICE-turn case, the improved path- following capability of the MLS 
system enables the M & S system to make better use of the fine-tuning potential 
of DICE maneuvers. 

As would be expected, there was no significant difference among the 
approach navigation systems studied when the traffic was all straight- in. Also, 
there was no significant difference between the ±60^ and ±40^ MLS systems for 
either the all-straight-in or the all-DICE-turn traffic sample. 

6.2.4 Summary of final-approach navigation canparison .- In general, for the 
particular ground-based, fixed-path speed-control M & S system evaluated, if 
the route structure contains a mixture of DICE-turn routes and straight-in 
routes, then MLS offers only a slight improvement over VOR/DME and ILS in the 
interarrival-error standard deviation performance. If the routes all contain 
two DICE turns, then the VOR/DME and MLS improvement over VOR/DME and ILS 
is more substantial, taking 2.6 sec off of the already fairly low 8.5-sec 
inter arrival-error standard deviation for the combined VOR/DME and ILS M & S 
system. For a geometry with only two-DICE-turn final-approach routes, the 
improved path- following capability of the MLS system over that of a conventional 
VOR/DME and ILS system enables the M & S control system to utilize some of the 
fine-tuning potential of the two DICE turns. For the M & S system studied, we 
can also say that wider coverage ±60*^ MLS system appears to offer no time- 
delivery precision over the basic ±40^ MLS system regardless of the arrival- 
route geometry. 


6.3 Wind-Error Effects 

With the ±60*^ MLS final-approach navigation system as a basis for compari- 
son, runs were made using three wind conditions to evaluate wind-error influence 
on this particular fixed-path M & S system. The three steady-state wind condi- 
tions were: (1) 10-percent error in wind strength and 5 percent error in direc- 

tion; (2) no error in strength or direction; and (3) -10 percent error in wind 
strength and -5 percent error in direction. The reference steady-state wind 
at ground level was 7.92 knots at 277^ (true north) with a linear increase in 
strength of 2.37 knots per 305 m (1 000 ft). Figure 18 shows that for the range 
of wind errors considered, the interarrival error standard deviation is not 
significantly changed. The dynamic schedule adjustment, which works to elimi- 
nate interarrival spacing violations, is responsible for keeping wind effects 
from building up. 


6.4 Tentative-Schedule Point 

The metering time used to determine if en route delay will be applied to 
regulate the traffic flow is actually not used by the terminal scheduler of the 
M & S system studied. Thus the metering is effectively uncoupled from the 
terminal sequencing and spacing. The ±60^ MLS system was used as a basis of 
comparison to determine the effect of moving the terminal tentative-schedule 
time from its normal 5 min to arrival fix. There appears to be no discernable 
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Figure 18.- Inter arrival delivery precision for various wind-error conditions. 
Final-approach navigation with +60° MLS; confidence interval of 95 percent 
for va 1 ues s hown . 
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Figure 19.- Inter arrival delivery precision for various tentative- 
schedule times from arrival fix. Final-approach navigation with 
±60° MLSj no wind error; confidence interval of 95 percent for 
values shown. 
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effect on the delivery precision, as shown by figure 19. This is because of 
the resequencing capability and the dynamic schedule adjustment which are trying 
to eliminate interarrival spacing violations. 

Although final delivery performance is not affected, the impact on the 
average flight time from the perimeter of the terminal to the outer marker is 
significant. Figure 20 shows that the average time of flight from the terminal 
perimeter is reduced frcxn 29.1 min to 19-4 min by moving the tentative sequenc- 
ing time from 5 min from the projected arrival (metering) fix time to 19 min. 

As figure 21 shows, the arrival-fix holding-stack delay time decreases as the 
time interval between the terminal tentative-schedule time and the metering fix 
is increased. Figures 22 and 23 illustrate this effect on the holding stack 
for identical traffic samples. 

Analysis indicates the average aircraft flight durations frcxn the 20-min 
ETA at the metering fix to runway touchdown are all about equal for the various 
tentative-schedule point runs made. Clearly, the closer the tentative-schedule 
time is made to the en route metering time, the closer we approach a coupled 
metering and scheduling system. Therefore, we can conclude that more of the 
delay for a given traffic arrival rate is taken at the more fuel-efficient 
higher en route altitude and less is taken at the arrival-fix holding stack if 
the metering and scheduling are coupled when both are time-based systems. This 
coupling should definitely be incorporated in the design of a metering and spac- 
ing syston using fuel-efficient profile descent paths. 

In addition to saving fuel, further safety benefits are gained by reducing 
the average terminal flight time. If delays are taken en route, the result 
is reduced terminal area traffic congestion and conflict risk. Since the 
terminal is the flight area of heaviest pilot work load, any reduction in con- 
troller directed maneuvers to delay and sequence aircraft is welcomed in the 
cockpit. Reducing the traffic volume within the terminal allows more controller 
attention time per aircraft, thus reducing controller stress. Another consider- 
ation is that higher en route rather than terminal delay reduces the interaction 
of high performance traffic with lower flying general-aviation visual flight 
rules (VFR) traffic. 


7.0 REAL-TIME STUDIES 

The primary purpose of conducting real-time pilot-in-the-loop experiments 
was to verify and demonstrate the fast- time results under representative cockpit 
work-load conditions. A simultaneous joint experiment was conducted to also 
study the effects of flight-deck systems operating in the M & S environment 
(ref. 1 3) . This situation subjected the pilot crew to a full range of cockpit 
tasks from cruise descent to runway touchdown. This prevented the crew, for 
example, fron concentrating on meeting delivery times at the outer marker at 
the expense of runway alignment and landing preparation. 
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Figure 20.- Average flight time from terminal boundary for various 
tentative-schedule times from arrival fix. Final-approach 
navigation with ±60^ MLS; no wind error; 35 aircraft per hour; 
confidence interval of 95 percent for values shown. 
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Figure 21 Average holding-stack delay time for various tentative-schedule 
times frctn arrival fix. Final-approach navigation with ±60° MLS; no 
wind error; 35 aircraft per hour; confidence interval of 95 percent for 
values shown. 
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Figure 22.- Traffic flow for 35 aircraft per hour data sample with 
a 5-min-to-arr ival-fix tentative-schedule time. Final- approach 
navigation with ±60^ MLS. 
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7.1 Ccx^kpit and Crew Considerations 

The TC7 Aft Plight Deck Simulator was integrated with the Terminal Area 
Air Traffic Model (TAATM) as described in section 3.5. A person acting as an 
air traffic controller working together with TAATM (see section 3.6) provided 
an environment and ATC commands to the simulated cockpit. The cockpit was 
treated by TAATM in the same way as software-generated aircraft controlled in 
the TAATM airspace. The task assigned to the research pilots was to fly an 
assigned STAR (described in section 3.6.2) into the simulated Denver terminal 
area with the onboard RNAV system. The control-system configuration used was 
a manual mode made up of the velocity control-wheel steering (VCWS) system and 
the autothrottle system with the pilot required to fly and land the aircraft 
manually. The cockpit is equipped with a set of electronic attitude director 
indicator (EADI) and electronic horizontal situation indicator (EHSI) displays 
(ref. 9) at both the pilot and copilot positions. 

For each of the experimental conditions studied, the equivalent of three 
complete two-man professional pilot crews were used. Three NASA research pilots 
were rotated as captain and first officer. These pilots are the crew members 
normally used for TCV flight experiments and therefore have considerable flight 
and simulator experience with the TCV aft flight deck system. 


7.2 Navigation and Route Factors 

As in the fast-time studies, the delivery performance of the ground-based, 
fixed-path speed-control system using VOR/DME and MLS was compared with the 
performance of this same system using current VOR/DME and ILS. Both ±40^ MLS 
and ±60^ MLS coverage were tested in real time to verify fast-time results on 
the issue of whether wider coverage MLS provided any benefit. As noted in sec- 
tion 2.0, the outer marker was moved 2 n.mi. closer to the runway for the MLS 
runs. In theory, the MLS's navigation precision should allow close-in runway- 
centerline intercept thereby reducing the final-approach ccmmon path. Later 
discussion of experimental results will indicate that this gate relocation may 
lead to real-world problems not originally anticipated. 

For the real-time case, time-of-ar rival errors at the outer marker for the 
TCV Aft Flight Deck Simulator can be independently and directly measured but 
the interarrival-error measurement must use a TAATM (real-time version) gen- 
erated aircraft as the "other" aircraft of a pair. Since verification is the 
goal, it would be desirable to isolate the real-time TCV simulator results from 
any TAATM fast-time effects. In addition, the aircraft inter arrival-error dis- 
tributions are not necessarily equal for the straight- in and DICE-turn routes 
of the geometry shown in figure 1 . With this situation there is a route inter- 
action effect in the interar rival-error measurement which, for the real-time 
single-aircraft case, cannot be directly measured in the same manner as was done 
in the all-TAATM, multiple-aircraft fast- time data runs. 

Appendix A details a procedure to obtain a strictly real-time TCV simulator 
estimate of the interarrival system performance. The procedure uses the piloted 
cockpit arrival-time-error variances at the outer marker for each type of route 
to calculate an equivalent-system interarrival delivery value. In addition to 
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the real-time arrival-error data, the procedure required knowledge of the rela- 
tive route loadings, which for this study were 40 percent for the straight-in 
route and 60 percent for the DICE-turn route. This approach yields the inter- 
arrival delivery variance or standard deviation which would be expected if the 
entire traffic sample was comprised of aircraft having the same delivery-error 
distribution as that obtained from the piloted-cockpit real-time data runs. 
Consequently, direct comparison with the fast-time inter arrival-error delivery 
precision is possible. 

Table 8 gives the number of initial real-time TCV-simulator data runs made 
for each approach navigation system and type of approach route. Each of the 
three crews repeated a particular experimental condition three times for a total 
of nine real-time data runs for that situation. 


TABIiE 8.- INITIAL REAL-TIME TCV AFT FLIGHT DECK SIMULATOR DATA RUNS 


Navigation system 

Number of runs for - 

DICE-turn 

routes 

Straight-in 

routes 

VOR/DME and ILS 

9 

9 

VOR/DME and ±40° MLS 

9 

9 

VOR/DME and +60° MLS 

9 

9 

1 


7.3 Real-Time Results 

7.3.1 Initial data runs .- Figure 24 gives a comparison of the fast-time 
and real-time interarr ival-delivery-precision results. Each of the first three 
real-time cases shown represents a total of 1 8 real-time experimental runs, 9 
of which were straight- in approaches and the other 9 were DICE-turn approaches. 
When the real-time results of the straight-in and DICE-turn runs are canbined 
using equation (A20) of appendix A, one unfortunately gets 8 degrees of freedom 
for the resultant estimate of inter arrival-error standard deviation. This pro- 
duces a fairly large confidence interval and less discrimination between output 
results from the F-test than if the route results could be directly pooled with 
16 degrees of freedom. Since there were no statistically significant differ- 
ences between the fast- time ±40^ MLS and ±60^ MLS systems or between those 
systems in real time, the real-time ±40^ mlS and ±60® MLS were combined to get 
a larger sample to represent the MLS system. The ±40® and ±60® MLS results were 
pooled at the route level and combined according to the procedure derived in 
appendix A. This combined MLS real-time result is shown on the far right of 
figure 24. Table 9 is a summarized comparison of the real-time and fast-time 
results. 
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Interarrival -error standard deviation, sec 



VOR/DME ±40° ±60° VOR/DME ±40° ±60° Both 

Its MLS MLS ILS MLS MLS and ±60 

Fast-time results Real-time results 


Figure 24.- Fast-time and real-time comparisons of interarrival 
delivery precision for various final-approach navigation 
systems. Confidence interval of 95 percent for fast- time 
data shown and of between 90.25 percent and 99.88 percent 
for real-time data. 
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TABLE 9.- SUB1MARIZED COMPARISON OF FAST-TIME AND REAL-TIME RESULTS 


Navigation system 

Standard 
of inter arrival 

deviation 
time error for - 

Real-time runs, 
sec 

Fast-time runs, 
sec 

VOR/DME and ILS 

6.3 

9.2 

VOR/DME and ±60° MLS 

1 2.3 

8.0 

VOR/DME and ±40° MLS 

1 5.7 

7.9 

Combined ±40° with ±60° MLS and VOR/DME 

1 3.6 

— 


The equal-tails F-test was used to determine whether there was discernible 
difference between the real-time and fast-time interarrival time-error variances 
at the 5 percent level. The hypothesis tested is 

^o* ® (fast time) “ ® (real time) 

against 

^1 • (fast time) (real time) 


where Hq is the null hypothesis; H] is the alternate hypothesis, and s^ is 
the variance. The number of degrees of freedom used were 400 for the fast- time 
data, 8 for the real-time VOR/DME and ILS , and 16 for the real-time combined 
±40^ and ±60^ MLS. A comparison of the fast-time and real-time results follows 
with equality indicated when the null hypothesis is not rejected. 


Fast-time VOR/DME and ILS = Real-time VOR/DME and ILS 
Fast-time ±40^ MLS ^ Real-time combined ±40^ and ±60^ MLS 
Fast-time ±60^ MLS / Real-time combined ±40^ and ±60^ MLS 


The real-time VOR/DME and ILS data seem to meet the initial objective of 
verifying the fast-time results. However, a significant factor is that not only 
do the real-time MLS results not agree with the fast-time data for the F-test 
but they also show a worse delivery performance. A comparison of the real-time 
data with the F-test indicates that even the VOR/DME and ILS real-time inter- 
arrival time-error variance was less than that observed for the real-time com- 
bined ±40^ and ±60^ MLS data. This real-time data comparison indicates there 
is only a 0.032 probability that the inter arrival time-error variance of the 
real-time VOR/DME and ILS is equal to that of the combined +40^ and ±60^ MLS. 


70 




in 


7.3.2 Gate location .- The question then is what factor is influencing the 
pilots* actions and overshadowing the accuracy of the MLS system? A check on 
the sequence of the real-time data runs eliminated a learning curve as the 
cause. The pilots themselves said they were not aware of a concious reaction 
to either the MLS or VOR/DME and ILS systems which would explain the real-time 
data discrepancy. A possible cause was the fact that the gate- to-r unway- 
threshold distance for the MLS system was less than that used in the VOR/DME 
and ILS runs. 

A look at the paths flown for the DICE-turn approaches indicated the 
following general tendencies: (1) The runway extended path seems to have been 
intercepted somewhat closer to the gate for the VOR/DME and ILS than for the 
MLS runs; (2) the turns on the MLS in general are somewhat more shallow with 
less bank angle than that of the VOR/DME and ILS; and (3) VOR/DME and ILS- 
system turns appear to more closely match the M & S algorithm assumed radius 
of turn. Figures 25 and 26 are comparable approaches flown by the same pilot 
illustrating these tendencies. Another set of 9 DICE-turn runs was made to 
verify these observations. These runs used the ±60^ MLS system but had the 
gate placed at the normal ILS positionr 12.04 km (6.50 n.mi.) from the thresh- 
old, rather than the 8.33 km (4.50 n.mi.) used in the regular MLS runs. The 
results of previous ILS straight-in runs were combined with the new ±60^ MLS 
(ILS gate) DICE-turn approaches to obtain an equivalent system-delivery perfor- 
mance. The interarrival-error standard deviation obtained was 9.7 sec and is 
plotted in figure 27 along with previous real-time runs for comparison. 

By using the P-test at the 5-percent level of significance, the hypotheses 

tested for the real-time results were: 

Hq* s (VOR/DME and ILS) = s (+60 MLS with ILS gate) 

against 

‘ ^ (VOR/DME and ILS) ^ ^ (±60 MLS with ILS gate) 

and 

^o* (combined ±40^ and ±60^ MLS) " (±60^ MLS with ILS gate) 

against 

^1 • (combined ±40^ and ±60^ MLS) ^ (±60^ MLS with ILS gate) 

Unfortunately, the null hypothesis could not be rejected in either case. Since 
the inter arrival-error standard deviation of the ±60^ MLS with its gate distance 
extended was between the standard deviation of the VOR/DME and ILS and that of 
the combined ±40^ and ±60^ MLS, the F-test indicated an equality to both systems. 
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Interarrival -error standard deviation, sec 


9 . 

k' 



VOR/DME Both ±40° ±60°MLS , 

Its and ±60° MLS ILS gate 


Figure 27.* Interarrival delivery precision of the real-time data 
runs for various final-approach navigation systems and gate 
locations. Confidence interval is between 90.25 percent and 
99.88 percent for values shown. 
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In searching further for confirmation of the possible gate- location effectr 
distances from the gate where the runway centerline was intercepted were measured 
from approach plots like those of figures 25 and 26 • The results of these mea- 
surements are summarized in table 10. The measured average intercept distance 
from the gate for the VOR/DME and ILS system was 0.39 km (0.21 n.mi.), whereas 
the average distance measured for the canbined ±40® and ±60® MLS was 0.89 km 
(0.48 n.mi.). The average values obtained seem to confirm the contention that 
the pilots intercepted the runway centerline extension closer to the expected 
M & S point (the gate) when the gate was farther from the threshold. However, 
instead of falling relatively close to one of the aforementioned values, the 
average distance measured for the new ±60® MLS run with ILS gate distance was 
0.67 km (0.36 n.mi.). This value is almost halfway between the values measured 
for the VOR/DME and ILS and for the combined ±40^ and ±60® MLS runs. At least 
these results are consistent with the interarrival-error standard deviation 
measurements because for those too the ±60® MLS run with ILS gate fell between 
the other two measurements (fig. 27) . 


TABLE 10.- MEASUREMENTS OF RUNWAY-CENTERLINE-INTERCEPT DISTANCE FROM THE GATE 


Navigation system 

Number 

of 

runs 

Average 

distance 

intercept 
from gate 

Standard 

deviation 

km 

n.mi. 

km 

n.mi. 

VOR/DME and ILS 

9 

0.39 

0.21 

1 .78 

0.96 

Conbined ±40° and ±60° MLS 

1 8 

.89 

.48 

1 .11 

.60 

±60° MLS with ILS gate 

9 

.67 

.36 

.80 

.43 


Although the magnitude of the average runway-centerline-intercept distance 
from the gate seems to vary with gate distance from the runway threshold, the 
small amount of real-time data does not give very good statistical discrimina- 
tion. For instance, the t-test for differences of mean used on the real-time 
results failed to indicate a definite difference in intercept distance. The 
hypothesis used was 


Hq! ^(VOR/DME and ILS) ~ ^(combined ±40° and ±60° MLS) 
against 

J ^(VOR/DME and ILS) ^ ’^(combined ±40° and ±60° MLS) 

where U is the average intercept distance from the gate. Because of the rela- 
tively poor resolution of the t-test resulting from the limited number of real- 
time runs, the null hypothesis could not be rejected even at a 1 0 percent level 
of significance. This same result is, of course, obtained when the average 
intercept distance for the ±60® MLS with ILS gate is tested against the inter- 
cept distance of the other two systems. Unfortunately hypothesis testing with 
the real-time data available does not permit a definite conclusion as to whether 
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the runway-center line- intercept distance from the gate is influenced by the 
gate distance from the threshold. 

Although the data is not as statistically conclusive as desirable, the 
distance of the gate from the threshold does appear to be a factor in pilot 
performance. The principal causes of outer-marker-delivery time errors are 
(1) A pilot does not start a DICE turn at the expected time; (2) the actual turn 
is not made with the bank angle assumed in the ground algorithm; and (3) the 
intercept, or merge, point with the extended runway centerline is at a different 
point than the gate location which was expected with the ground ETA predictive 
algorithm. The latter two situations appear to be the contributing factors when 
the gate is moved closer to the runway. A plausible explanation for the worst 
delivery performance experienced when the gate was moved closer to the threshold 
is that the primary concern of the pilot at sane point in his approach becomes 
that of configuring and aligning the aircraft for executing the normal, safe 
IFR landing required in the real-time simulation runs. Given that the pilots 
had experience with the EHSI map in normal TCV landings, the subconscious reac- 
tion may be to align and configure their aircraft for landing at sane distance 
fron the threshold which training and experience tells them is required for an 
acceptable touchdown. Another consideration is that for a given glide slope, 
runway-alignment maneuvers must be performed closer to the ground when the gate 
is moved in. All of this considered, we submit that in this investigation, the 
MLS gate was moved close to or inside the point where the inherent pilot concern 
for landing impacted the M & S system capability for precise time delivery. 

Although the gate-location effect is not predictable in the fast-time study, 
this effect is an interesting and significant finding in view of the fact that 
the next generation of commercial aircraft will be equipped with EHSI displays. 
When the MLS system gets implemented, the possibility of moving in the gate to 
shorten the common IFR path will be considered (ref. 1 4) . Conbine these factors 
with the pressure to achieve maximum IFR capacity through precise time delivery 
and the same situation as that experienced in the real-time runs will exist. 
Although the experimental findings gave an indication of what might be expected, 
it is possible that further experience and training might modify the impact on 
delivery performance of moving the MLS gate closer to the threshold. 

The effect of an onboard coupled control system in removing pilot variation 
is a related question. The limited results obtained raised some interesting 
and important issues which should be addressed in further studies. 


8.0 SUMMARY OF RESULTS 

The fixed-path (speed-control) metering and spacing (M & S) system described 
in this report assumes an aircraft two-dimensional area navigation (RNAV) capa- 
bility. The chief objectives of this investigation were to evaluate the delivery 
performance of an earlier proposed M & S system and to determine if using the 
microwave landing system (MLS) for approach guidance improved the system's 
delivery performance over that of very high frequency omnidirectional range/ 
distance measuring equipnent and instrument landing system (VOR/DME and ILS) . 

The primary performance measure used was the aircraft inter arrival-error stan- 
dard deviation at the outer marker. Facilities used in the study are the 


75 


Terminal Area Air Traffic Model (TAATM) and the Terminal Configured Vehicle 
(TCV) Aft Flight Deck Simulator. Both fast-time (batch-mode) and real-time 
pilot- in- the-loop experiments were conducted. Results of the investigation are 
discussed in sections 5.0, 6.0, and 7.0. The following is a list of the major 
findings of the investigation: 

1 . The fast-time simulation data show that the M & S system aircraft 
inter arrival-error standard deviation at the outer marker is only reduced fran 
9.2 to 7.9 sec when the ±60^ MLS replaces the VOR/DME as the RNAV data source 
during final- approach, fine-tuning, DICE-directed maneuvers. Although this is 
true for the particular geometry studied, further investigation indicated this 
result reflects the interaction of two types of approach routes. For only 
straight-in routes there was no verified difference between the interarrival- 
error standard deviation of the MLS and ILS systems; however, there was a more 
substantial 2.6-sec advantage for the MLS system if the geometry contained only 
DICE-turn final- approach routes. The DICE-turn routes have the downwind-base 
configuration normally used at most terminals. 

2. There was no statistically significant difference between the time- 
delivery performance of the ±40^ MLS system and that of the ±60® MLS system. 

For the particular M & S system studied, the wider coverage ±60^ MLS appears 

to offer no time-delivery precision over the basic ±40^ MLS system. This would 
appear to be true regardless of the arrival-route geometry as long as the DICE 
turns are made within MLS coverage. 

3. Real-time VOR/DME and ILS data results were consistent with the fast- 

time data runs but the real-time MLS results indicated a worse delivery per- 
formance not only compared with the fast-time data but also compared with the 
real-time VOR/DME and ILS system as well. The gate location for the ILS runs 
was placed 12.04 km (6.50 n.mi.) from the threshold. Since the increased 
accuracy of the MLS system allows precise final intercept, the gate location 

for the MLS runs was moved 3.70 km (2.00 n.mi.) closer to the threshold in order 

to reduce the final-approach common path. Though not as statistically conclu- 
sive as desirable, the real-time data indicate that moving the gate 2.00 n.mi. 
closer to the runway threshold may affect the pilot's implementation of final 
DICE-turn maneuvering needed to achieve M & S delivery precision. It appears 
that if the gate is moved significantly closer to the threshold than normal, a 
pilot with the displays available in the TCV aft- flight-deck cockpit may tend 

to configure and align the aircraft for landing at the same distance he normally 

does rather than precisely following the M & S-derived fine-tuning DICE instruc- 

tions. This raises some important issues concerning the interactions of 
electronic horizontal situation indicator (EHSI) displays, close-in final inter- 
cepts, and M & S systems which need further investigation. 

4. The en route metering function of the M & S system studied, even though 
accanplished with a time-based process, only regulates the aircraft arrival rate 
and operates independently of the terminal schedule process. That is, the ter- 
minal scheduler has no knowledge of the metering-scheduled outer-marker time and 
establishes its cwn scheduled outer-marker time when the aircraft reaches the 
terminal's tentative-schedule point. This study demonstrated that most of the 
arrival (metering) fix holding is transferred to more fuel-efficient, higher 
altitude en route delay as the terminal's tentative-schedule point is moved 
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closer to the meter ing -process time point. Clearly the closer the terminal 
tentative-schedule time is moved to the en route metering time, the closer a 
coupled metering and terminal scheduling system is approached. Therefore, 
coupling should be incorporated in the design of a metering and spacing system 
utilizing relatively fixed, fuel-efficient, profile descent paths. Also, a 
reduction in terminal delay not only saves fuel but increases safety by reducing 
congestion and lowering the work load of both pilots and controllers. 

5. The M & S system used in the investigation had several limitations and 
constraints which are discussed in section 5.0. For example, the ground rules 
used to assign holds tended to be conservative and inflexible, thereby creating 
traffic-flow gaps unnecessary for safety. Another constraint noted is the lack 
of information in the scheduling algorithm on the path-shortening and stretching 
capability available. To achieve a relatively fixed-path configuration this 
information, by design, is used only by the schedule-maintenance algorithm as 
a backup to realize a scheduled marker time. This procedure places highest pri- 
ority on achieving the scheduled marker time by speed control with secondary 
priority on holding aircraft- and resequencing depending on the aircraft loca- 
tion. The limitations of this process are compounded by the holding limitations 
mentioned previously. If the scheduler was given more knowledge of the control- 
lability possible with path shortening and stretching, then speed control could 
still be given highest priority but the secondary priorities would have path 
change over holding. This should yield smoother traffic flow but at the expense 
of some variability in aircraft approach paths. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 6, 1981 
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APPENDIX A 


SYSTEM INTERARRIVAL TIME-TERROR PERFORMANCE FROM 
SINGLE-AIRCRAFT, ROUTE -DEPENDENT ARRIVAL ERRORS 

During the real-time experiments, the TCV Aft Flight Deck Simulator was 
treated by the Terminal Area Air Traffic Model (TAATM) as merely another air- 
craft in its airspace to be controlled. The time-of-arrival error at the 
outer marker for the TCV simulator can be independently and directly measured 
but any interarrival-error measurements must use a TAATM (real-time version) 
generated software aircraft as the "other" aircraft of a pair. Since one of 
the objectives of the real-time investigation was to verify fast-time TAATM 
performance, it would be preferable to have a "pure" real-time result not 
affected by the quantity to be verified. In addition, the arrival errors are 
different for the straight- in and DICE-turn routes used in the geometry. There- 
fore, there is a route interaction in the interarrival-error measurements which 
cannot be directly measured in the single-aircraft, real-time data case as was 
done in the multiple-aircraft, fast-time data runs. The development in this 
appendix allows the measured, piloted-oockpit error variance at the outer marker 
for each type of route to be converted to an equivalent-system interarrival- 
error variance. This interarrival-error variance is what would be expected if 
all aircraft in the terminal system had the same time-error characteristics as 
those measured for the TCV simulator. 


Arrival Error of Single Aircraft 

This section will show how to determine the single-aircraft arrival error 
when the arrival could be from either of two types of routes, each with its own 
delivery-error distribution. The arrival errors for each of the types of routes 
are independent of each other. 

Let the route be defined by variable R so that R*| has an aircraft time 
error of delivery distribution x such that f(x|R-|) - N(u-|,a-|2), the notation 
being that the conditional distribution of x given route 1 is a normal distri- 
bution with mean u-| and variance a-|2. Similarly, R 2 has an aircraft time- 

error distribution x such that f(x|R 2 ) ~ N(U2/C^2^) with probability P of 
each route as follows: 


P(Rl ) = Pi (A1 ) 

P(R2) = P2 (A2) 


and 


Pi + P 2 = 1 (A3) 
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The following sketch illustrates the conditional time-error distributions for 
each type of arrival route: 



We can now state the following 


R 2 

fx(x) = 53 fx(x|Ri)P(Ri) 

i=1 


f X ^ “ P] fx ^ ^2 I ^2 ^ 


fx (x) 



(x-ui ) 2 


(X-U2) ^ 

1 

2ai2 

1 

202^ 

Pi 1 — e 

Oi ^ 


+ P2 ;= e 

a2pn 



(A4) 
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The mcment-generating function for x, represented by M^Ct), is determined 
by the following: 


■i: 


Mx(t) = E{e^^} = I f(x) dx 


^ tu.r ■' t(x-u ) 

= P*l e II := e ' e 


(x-Ut )2 


2a I 


ai \]^ 


dx 


tu« ^ t (x-u ) 

+ P 26 2 I ;= e ^ e 


/ 

_o 


(x-U2)2 


2a-] 


O'] i/ItT 


-■ dx 


(A5) 


where E{e*-^} is the expectation of e^^. By using the known results for 
normal distributions we get 


M»(t) . + p^^ltU2-KG2V)/2l 


(A6) 


By using the property of M^Ct), 


E{x) = 


= PlUi + P 2 U 2 


(A7) 


t=0 


and 


e{x 2} = [Mx(tn 

dt2 


=, P] (a-i 2 + ui 2) + P2(a22 + U2^) 


t=0 


(A8) 
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II 

I 

We can now obtain the variance of x (i.e./ Var{x}): 

Var{x} = E(x2) - [e(x)]2 

= Pt (a-^2 + ut^) + P2(a2^ + U2^) - (PiUi + P 2 U 2 ) ^ 

Var{x} = Piai2 + P 2 Q 22 + P] P 2 (Ui - U 2 ) ^ (A9) 

Interarrival Error of Two Aircraft With Known 
Arrival-Error Distribution 

Given there are two particular aircraft arrivals such that 



where the first aircraft A has target time t^ and error x^ and the second 
aircraft B has target time tg and error X 3 . 

The errors x;^ and xg are independent: 


^A ^ ~ ^x I 

Xb ~ N(UB,aB2) = fx(x|B) 


(AT 0) 


Therefore, the interarrival error y between A and B is 


y = [(tg + xb) - (t^ + x^)] - (tB - t^) = Xb - xa 


(All ) 


We can thus say the following: 


E{y} = Xg - Xa = Ub - Ua 


(A1 2) 
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and 


Var{y} = E{y^} - [E{y}]2 

= e{ (xb - - [e{xb - *A^^^ 

= e{xb^} - Ub^ + e{xa^} - 
= Var{xB) + Var{x;^} 

Var{y} = Op2 + (AT 3) 

The interarrival error y, given the order and arrival error of each aircraft 
is 


fy(y|A,B) ~ N(Ub - Ua, 0^2 + 0b2) 


(A14) 


General Interarrival Error of Two Aircraft 

We can now determine the interarrival time error between two aircraft when 
either aircraft could come from either of two types of routes with each route 
having its own error distribution. With the symbols already defined and from 
the following sketch, a new randcxn variable C is introduced such that 


C = 1 
C = 2 
C = 3 
C = 4 


(A = Ri , B = Rt ) 
(A = Ri , B = R2) 
(A = R2, B = Ri ) 
(A = R2, B = R2) 
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For independent arrivals the probability P is defined by; 
P (C = 1 ) = P] 2 ^ 

P(C = 2) = PiP2 
P(C = 3) = P]P2 ^ 

P(C = 4) = P2^ 


If the inter arrival time error between two aircraft A and 
y, then 


fy(y|c=i) ~ 


fy(ylc=2) ~ 


fy(y|c=3) ~ 
fy(y|c=4) ~ 


The moment-generating 
the following: 



N(0, 2ai2) 


f 2 ^ 

N(^Vic,2r ^c,2; 

1 ~ N(U2 - Ut , Ot 2 + 

022) 

N^Uc,3, 

~ N(Ui - U2, Ot 2 + 

022) 

n(o, ~ 

N(0, 2022) 


function for 

y^ represented by 

My(t) 
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20 


(A15) 


is defined by 


(A16) 


is determined by 



dy 

(AT 7) 
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By using the known results from normal distributions we get 


My(t) = Pi 2 e 


t^c,2+(c^c,2t^)/2 

‘ '' + P] P 26 •- ^ '' -i 


+ P2Pie 




+ P2^e 


,a-,2t2 rtU2-tUi+(aT2+a22)t2/2l 

My(t) = Pi2e 1 + PTP2e*- 


+ P2Pie 


[tUT-tU 2 +(a-, 2 + 032 ^ t 2 / 2 j 


2.2 


P2^e 


2^°2 ^ 


(A18) 


(AT 9) 


By using these properties of My(t), 


E{y) - - 0 


and 


E{y2} = ^[m (t^l 

dfc 2 Ly Jlt=C 


we can say the following: 


E{y2} = Pi 2(2ai 2) + 2PiP2 




(U2 - Uij 


+ P2^{202^] 


(A20) 


Since E{y} = 0, 


Var{y} = E{y2} 


(A21) 


Confidence Interval of Interarrival-Error Standard Deviation 

This section shows the determination of some interval around the sample 
estimate of interarrival-error standard deviation 0y defined by end points 
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(confidence limits) for which sane measure of confidence can be assigned. The 
terms ^H,l lower and upper bounds of the confidence inter- 

val for a-| determined fran the chi-square such that 



^L,l 

= t>L,lcri 


C^H,1 

= 

and 

and 

*^H,2 3' 



= bL,2C^2 


^H,2 

= bn,2<^2 

where t>L 

and 

b^ are 


degrees of freedom. Therefore, 


(A22) 


(A23) 




C'l,! = <^1 = 

and and and 
^ L ,2 = ^2 = ^ H ,2 

I 1 I 

°L,y = = °H,y 


Arrival-error standard deviation 
confidence bounds 

Transform function (eq. (A20) ) 

System inter arrival-error confidence bound 


with cJl,y ^H,y defining the confidence interval around Since all 

terms are positive, the transform function for the lower bound of Oy is 
obtained by substituting the lower bound for each of the routes (i.e., ' 

Ql, 2 ) into equation (A20) to get 


<^L,y = {Pl2[2(bLjaT)2] + 2 Pt P2[(bL,iaT ) 2 + (bL,2^2)^ + (^2 - U] )C1 

+ P22[2(bL,2«^2)^]>^^^ (A24) 


A similar substitution yields the upper bound. If the degrees of freedom are 
the same for the data fron both routes, then 


t»L,l = ^l,2 * 


(A25) 
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By substituting into equation (A24), 


^L,y 


= br 


;2Pi 2ai 2 


2P^P2 





+ 


2P2^a2 



Rewriting terms of Oy fron equation (A20) yields 


(A26) 


'L,y 


= br 


0 2 _ 
<Jy 


2PtP 2(U2 - Ui ) 



(A27) 


A similar process for the upper bound gives 


*^H,y - 



2PtP2(U2 - Ut ) 




(A28) 


A procedure to define the end points of a confidence interval has been 
determined. Now some measure of confidence is needed for the interval. The 
following sketch illustrates the confidence limits: 



If P = 0.9500 that G-j is within the interval ^ Gj ^ g^ j and also the 

same probability that 02 is within the interval = ^2 = then 

P = 0.9025 that Gj and C 2 are both within their respective confidence bounds. 

If both ^ ai ^ Ohj and = <^2 = <=^H,2 are true, then ^ Oy ^ Qg „ 

is also ture. See eq. (A24).) If 
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P = 0.025 that 0*1 > 0 h,T same probability that 0] < 0 l^*| 


and 


P = 0.025 that 02 > ^h, 2 same probability that 02 < 


then 


P = 6.25 X 10*"^ that both 0] > 0h, 1 ^2 ^ ^H,2 same 

probability that both o-j < 0 lj and 02 < 0l,2 

Therefore/ P = 1.25 x 10“^ that a-j and 02 are both above or both below 
their respective confidence bounds. From equation (A24)/ when these conditions 
exist 0y is outside its confidence bounds defined by ^H^y ^L,y 

result of this exercise is a bounded probability interval for the confidence 
interval of Qy. If (o^L/ 1 '^H/1 ) (^L/2'^H/2) <3efine 95.000-percent con- 

fidence intervals for o-j and 02 / respectively/ then 0y will fall between 
^L,y *^H/y somewhere between 90.250 to 99.875 percent of the time. This 

gives a feeling for the confidence which can be placed on the interval defined 
by c?L,Y «^H,y 


Fast-Time Experimental Verification 

The two types of routes used in the system studied were the straight-in 
BYERS route and the two-DICE-turn approach of all the other arrival routes. 

As a verification for equation (A20) , the aircraft arrival-error standard 
deviation was obtained from a fast-time/ no-BYERS-tr af f ic data run and also for 
an all-BYERS sample. These arrival-error standard deviations were combined in 
equation (A20) to obtain the calculated system interarrival-error standard 
deviation. Figure Al shows this calculated value was indeed quite close to the 
measured system interarrival-error standard deviation of the normal-system/ all- 
route data run. 
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data run all -strai ght-in and 

al 1 -DI CE-turn traffic run 


Figure Al Interarrival delivery precision frcwi a measured, normal 
all-route data run and calculated value from arrival precision 
of an all straight-in and all-DICE-turn traffic sample. Final 
approach with ±60^ MLS navigation; no wind error; confidence 
interval of 95.00 percent for measured values shown and between 
90.25 and 99.88 percent for calculated value. 
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SPEED-CONTROL ALGORITHM (SPDOPT) 

The nominal speed on a given flight-path segment for an aircraft is the 
highest acceptable speed for that aircraft dependent on its position within the 
terminal area and its performance class. The schedule-maintenance logic of the 
M & S system, however, is capable of selecting an alternate speed if a speed 
option has been specified for the segment and if the control logic has deter- 
mined that an alternate speed should be chosen to achieve a desired SMT. Assum- 
ing that Vj^ and are the permissible high and low indicated airspeeds 

for a given segment, the speed-control algorithm is responsible for choosing 
an indicated airspeed Vjg such that ^ Vjg « 

This appendix details the equations used to calculate a ground speed Vs^ 

equivalent to Vj when speed control is applied on a segment for a given air- 
"*■0 

craft. The velocity limitations on the segment and all succeeding segments on 
the route are presented to the algorithm as indicated airspeeds. In order to 
perform estimated-time-of-arr ival calculations, however, the algorithm first 
converts these velocities to equivalent ground speeds for the altitudes and 
wind speeds which are predicted to be encountered by the aircraft. The desired 
ground speed Vg^ calculated with the algorithm is also converted to its equiv- 
alent indicated airspeed before delivery as a speed command to the aircraft. 

The following discussion includes the equations used to convert indicated air- 
speeds and accelerations to their equivalent ground speeds and accelerations 
as well as a derivation of the equations used to calculate the desired ground 
speed Vg^ . 


Relationship of Indicated and Ground Speeds/Accelerations 

By using a second-order approximation, ground speed Vg is defined by 
the following; 


where 

Vg = k(l + ah + Bh2)Vj - w 

h 

altitude 

Vi 

indicated airspeed 

W 

wind vector 

k 

crab-angle an descent-angle correction factor 

ot,3 

constants 
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The approximate ground acceleration a-] for segment 1 S-| is 

a"! = aj(l + ctho + + bVj^ (cx + 2BhQ) 

where/ for k ^ ^ and 

dVg 3Vg dh 9Vg dVj 

dt 3h dt 3Vj dt 


the following definitions apply; 


(B2) 


aj constant indicated-airspeed acceleration (negative acceleration for 

dVj 

arrival case) , 


dh 

b constant descent rate, — 

dt 


a = 1 .48060 X 1 0“5 


6 = 9.071 43 X 1 0-'''' 


hQ initial altitude at the start of S-| 

initial indicated airspeed 


Likewise, where h-) is the initial altitude at the start of segment 2 
S2, the approximate ground acceleration 32 for S2 is 


32 — (1 ®hi + Bh] + bVj^ (ex + 2Bhi ) 


(B3) 


where Vj is the maximum possible indicated airspeed on S-| . Similarly, where 

h2 is the initial altitude at the start of segment 3 S3, the approximate 
ground acceleration 33 for S3 is 


33 = ai(T + *^h2 + Bh2^) + bVj 2 ((X + 2Bh2) 


(B4) 


where Vj2 is the maximum possible indicated airspeed on S2* 
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Definition of Terms Used in No-Speed-Overlap Case on Succeeding Segments 
Figure B1 illustrates the no-speed-overlap case for segments Sjj for 


n = 1 , 

and 3 where 


Vil 

initial ground speed on S-j 


Vhi 

highest possible ground speed on 

Si 

^1.1 

lowest possible ground speed on 

Si 

Vh2 

highest possible ground speed on 

S2 

CM 

> 

lowest possible ground speed on 

S2 

Vh3 

highest possible ground speed on 

S3 


desired ground speed on S] 




Distance 

Figure B1 .- No- speed-over lap case on succeeding segments. 


For this case ^ ^ ^ ^ Vh 2 ^ cases ^ or 

Vjj^ is set equal to . The determination of Vg^ requires the expression 

of XLLT(Vg^ ) and ETAR(Vg^ ) as explicit functions of Vg^ . The following 
two sections of this appendix define XLLT(Vg^) and ETAR(Vg^ ) for this case. 
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Computation of XLLT(Vs^) for No Speed Overlap on Succeeding Segments 

The term XLLT(Vs, ) is defined as the time of arrival at the OM if S-j 
is flown at and the other segments are flown at their lowest respective 

velocities. Figure B2 illustrates the time, distance, and velocity relation- 
ships used to determine XLLT(Vg^). We can say 


XLLT(Vg^ ) = XLLT(Vji^ ) + AXL^Tg^ + AXLLTg2 


(B5) 


where 


AxLLTg^ = (ti + t2) - (t^ ' + t2’) 


and 


(B6) 


AxLLTs2 " ^^3 + t4) - (t^' + t^') 


(B7) 


U 

o 



Figure B2.- Time, distance, and velocity diagram to illustrate XLLT(Vg^ ) . 
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The value XLLT(Vjj^) is defined in a manner similar to XLLT(Vs^) with the 
exception that S-| is flown at Vpj^ . This value is calculated external to the 

speed-control algorithm and is presented as an input constant to the routine. 

By referring to figure B2 and assuming constant acceleration, we can write the 
following equations: 


to = 


Vht - Vil 


ai 


(B8) 


t] = 


Vst - Vht 


(B9) 


d2 (S] - dg) - di 

t2 = = 

Vst Vst 


(BIO) 


do = 


Vm + Vi^\ - Vi^2 


to = 


2 a") 


(BIT) 


di = 


Vst + VhA Vs^ 2 _ Vh^ 2 


It] = 


2a-| 


(B12) 


t * + t * = 
1 ^2 


d] + d2 St - do 


Vht Vht 


(B1 3) 


Therefore, 


Axllt. 


N - ''H,' 


Si 


V 


ai 


+ St - 


Vti 2 _ w . 2^ 

^Ht '^IT \l 1 1 


'Vst^ - Vht^' 


2ai /\Vs^ Vh^/ V 2aTVg, 


(B14) 


Again referring to figure B2 we can write 


Vl2 - Vst 


t3 = 


32 


(B15) 
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<34 

t4 = 

Vl2 


(B1 6) 


t ' + t ' = 
^3 4 


Vl2 - Vht 


32 


(B17) 


d3 = 


^Vl 2 + Vs^\ - Vs^ 2 


t3 = 


232 


(B18) 


<33 + d 4 =1 


^VL2 + 


— y(t3’ .t,., - 


Vl2^ - Vhi 2 


232 


(B19) 


d4 = 


VSt 2 - Vht 2 
232 


(B20) 


Therefore, 


AxllTs2 = 


Vht - Vsi 


32 


Vsi^ - Vhi^ 

232VL2 


(B21) 


Computstion of ETAR(Vg^ ) for No Speed Overlep on Succeeding Segments 

The term ETAR(Vg^ ) is the estimated time of srrival 3 t the outer msrker 
using Vg^ on S-| and the highest velocities permitted on the remaining seg- 
ments. Figure B3 illustrates the segment-speed relationship used to determine 
ETAR(Vg^ ) . We can write ETAR(Vg^ ) as 

ETAR(Vg^) = ETARCVg^) + AETARg^ + AETARg 2 (B22) 

The expressions for AETARg^ and AETARg 2 sre similar to those used for 
AXLLTg^ and AXLLTg 2 with V ^2 replaced by Vjj 2 * By using equation (B1 4) 
and by replacing equation (B2lf, then 
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Solution of Speed Control for No-Speed-Overlap Case on Succeeding Segments 
The criteria used to determine the segment velocity Vg^ are the following 


and 


SMT = XLLT' (Vg^ ) 


(B25) 


SMT = ETAR(Vg^ ) + k' [FDELAY(Vg^ ) ] 


(B26) 


where k* is empirically defined as 


k* = k + (1 - k) 


/vm 

- VsA _ 

n 

If _i_ 

- k)vm 

4_ 

-0 - k) 

Wi 

- VlJ ■ 

i\ 1 

Vh^ 

-V.,J 

T 

1 

< 


= + k2Vg^ 


(B27) 


where k-| and k2 are the bracketed terms in the preceding equation. The 
value XLLT'(Vg^) is defined as the outer-marker arrival time if S-| is 
flcwn at Vg^ and the remaining segments are flown assuming a percentage of 
the future delay capability FDEIjAY(Vg^) will be used in order to achieve the 
desired SMT. This percentage approaches 100 percent as Vg^ approaches . 

The value FDELAY(Vg^) is determined with the following equation; 


FDELAY(Vg^) = XLLT(Vg^) - ETAR(Vg^ ) 


(B28) 


The following drawing illustrates these functions: 
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If equations (B5) , (B1 4) , and (B21 ) are applied to equation (B28) for XLLT(Vs^ ) 
and equations (B22), (B23) , and (B24) are used for ETAR(Vg^ ) , then 


FDELAY(Vg^ ) = 


VHi ^/1 1 ' 

PDELAY(Vh. ) 

' 2a2 \vl2 Vh2^ 


j_/_L _L' 




= PDt + FD2^Vs^2^ (B29) 

where FD-] and FD 2 are the bracketed terms in the preceding equation. If 
we substitute equations (B23) and (B24) into equation (B22) , the result is 


ETAR(Vq, ) = 

1 V- 


Si 


Vi 2 - 



+ 


Si 

ETAR(Vh, ) 


'Vht ^ - Vi 


2 _ V. 2 \ Vh^2 ^ ^ 

2aiVHi / 2a2Vu2 ^ 


+ VSi 


1 1 

2 ai a 2 


+ Vg 


2 a 2 VH-: 


ET] 


+ ET2 + ET3(Vg^) + ET4(Vg^2) 


(B30) 


where ETt , ET 2 f ET 3 , and ET 4 are the bracketed terms in the preceding 
equation. 
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Substituting the final forms of equations (B27) , (B29) , and (B30) into 
equations (B25) and (B26) yields the following: 

(Vg^ ) SMT = ETt + ET 2 (Vg^ ) + ET 3 (Vg^ 2) + ET 4 (Vg^ 3) + (Vg^ ) 

+ k2FDi(Vg^2) + k^FD2(Vg^3) + k2FD2(Vg^4) (B3T ) 

which can be written in the form of a polynomial in V3^ : 

AVg^ 4 + BVg^ 3 + CVg^ ^ + DVg^ + E = 0 (B32) 

where 


A = k2(FD2) 

B = ET4 + k-| (FD2) 

C = ET3 + k2(FD*i ) 

D = ET2 + k-| (FD^ ) - SMT 
E = ET] 

Solution of Speed Control for Speed-Overlap 
Case on Succeeding Segments 

Thus far in the discussion we have only considered the situaton for no 
speed overlap on adjacent segments, that is, ^ ^H2* overlap case 

in which the speed ranges are equal on segments 1 and 2, that is, 

and = ^L2' ^ restriction must be placed on as a function of 

in order to determine a desirable V3^ . The ground rules for the M & S system 
do not allow for increases in speed on succeeding segments. Therefore, whenever 
^S-| ^ ^H2' highest velocity permitted on segment 2 must be equal to V3^ , 

that is Vjj2 set equal to V3^. Figure B4 illustrates this speed overlap case 
for segments for n = 1 , 2, and 3. By going through a similar development 

as was done for the no-overlap case, a polynomial in V3^ can be determined: 

AVg^ 4 + BVg^ 3 + cvg^ 2 + ^Vg^ + E = 0 (B33) 
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Distance 


Figure B 4 .- Speed-overlap case on succeeding segments 

where 

A = k2(FD5) 

B = ki (PD6> + k2(FD5) 

C = ET7 + k-| (FD5) + k2(FD4> 

D = ETe + ki (FD4) + k2(FD3) - SMT 
E = ET5 + ki (FD3) 


for which the following apply; 
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1 + 0th2 + 3h2^ 

Cl = 

1 + ahi + 3hi 2 


■=1Vh,2 sj 

FD3 

233 Cl 


h Ci\ Vg^2 

FD4 = FDELAY(Vht) + Vhi 

' V^2 ^3/ 2a2VL2 



1 

32 


FDg 


1 

232VL2 


ET5 = Si 



+ — + 


Cl 2ai 


C]Vh^2 


2a 3 


ETe = ETAR(Vh^ ) 



S2 

Si + 

, 


1 

ET7 = — 
233 


c 


1 


233 


7 00 


S2 


^Hi^-Vil 

2^1 
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GENERAL DESCRIPTION OF RNAV 
Position and Track Information 

In order to obtain the distances and angles needed in the RNAV equations, 
the aircraft estimated position must first be determined- With defined 

as the actual aircraft position and (Xy,Yv) as the VOR/DME location, a line 
can be drawn from the VOR/DME station to the aircraft, as shown in figure Cl (a) - 
The angle 0 of this line from north N and its length p are determined 
from 


P = ^(Xa - Xv) 2 + (Yft - Yv) ^ 



If Ap and A6 are the ccxnbined airborne and ground errors, then the aircraft 
location errors (Ax,AY) are determined from 


Ax 


-p sin 6 cos 0 


ag' 

Ay. 


p cos 0 sin 0_ 


_AP. 


(C2) 


The location errors are added to the aircraft's position to simulate the raw- 
navigation measured location : 


Xm = Xa + 
Ym = Ya + 



(C3) 


The measured or input cross-track error D^j^ can now be determined. A straight 
line is drawn between the two waypoints defining the current aircraft path as 
shown in figure Cl (a) . The distance S from the measured aircraft position 
to waypoint 2 is 


S 


^(Xm - Xp,2)2 + (Ym - Yp,2)2 


(C4) 
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and, the bearing angle 6 is 


6 


7T 

2 


tan"^ 


/ Xp,2 - M 

V^P,2 - 


(C5) 


The bearing angle from waypoint 1 to waypoint 2 is determined frcxn 


'i'd 


ir 

2 


tan~'> 


/ Xp,2 ~ 

W,2 - Yp,l/ 


Angle r) is determined from 


n = 6 - 'Fd 


(C 6 ) 


(C7) 


The perpendicular offset of the aircraft's measured position from the course 
defined by waypoints 1 and 2 is the cross-track error and is given by 


= -S sin ri 


(C 8 ) 


If part of figure Cl (a) is slightly redrawn in figure Cl (b) with an aircraft 
ground velocity Vg and a ground-track angle ^ 3 , then a track-error angle 
is defined as follows: 




The cross-track deviation rate b is 


b = Vg sin 


(CIO) 


In order to smooth the raw navigation input data, a simple complementary- 
type filtering is done on the track position. If the following are defined as 

D aircraft position distance off track 

navigation-measured distance off track given by equation (C 8 ) 

estimated position of D at time n after smoothing 
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Dn track-rate information at time n 

predicted position at n using D^-i and 
then the smoothing function used is 

Dn = Dp,n + Kf (Din ~ Dp,n) = ^n-^ + + KftOin “ (Dn-1 + ^tbn)] (Cn ) 


where At is the position update time and Kf is a constant selected at 0.25. 
The filter action is shown in the following sketch: 


D. 


1 n 



Waypoint 1 


Waypoint 2 


The term defines the smoothed estimated position of D at time n, which 

is shown as (X 3 ,Y 3 ) in figure Cl (b) . 

Figure C2 contains the horizontal path control algorithm for RNAV guidance 
(see ref. 15), for which the following apply: 

desired bank angle 

true airspeed 

R radius of curvature when flying a curved path 
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’^NOM 


tZ0° when in RNAV turn 
0® otherwise 



7.12 


0.5 for IAS < 100 knots 

= 0.68 - 0.0018 IAS for 100 < IAS i 300 knots 

0.14 for IAS > 300 knots 


L 


V K. sin 

a y 


I 


1 30® for X < 30® 
90® for X > 90® 
X otherwise 


X 


470 Id 



30® 


Figure C2,- Horizontal path-cx>ntrol algorithm for RNAV guidance. 
(1 ft = 0.3048 m; constants will change if metric values are 
substituted. ) 
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Turn Anticipation 

For an RNAV aircraft flying from waypoint A to B to C as in fig- 
ure C3, the objective of turn anticipation is to start the aircraft turning 
toward waypoint C at a point F such that the aircraft flies a smooth path 
toward segment BC with no overshoot of the course. As the aircraft flies frcan 
A to B (fig. C3) the following calculation is made to determine the distance 
d from B to initiate roll-in to the nominal bank angle '^’nOM the turn: 


d = 



(Cl 2) 


where A'l' is the change in heading between segments AB and BC, R is the 
radius of curvature defined by 


R = 



g tan <I>nom 


(Cl 3) 


for “ 20°, and r^ is the roll-in distance to establish the nominal bank 

angle: 


ri = AtVa 


(Cl 4) 


with At as the time required to roll into nominal bank: 


"^NOM 

At = -r (Cl 5) 

^Limit 


where ^Limit “ ^ degrees per second. A roll-out distance is also established 
to make a smooth merge with course BC. A normal bank angle of 0^ is desired 
at the start of the straight segment if the aircraft is on course and track. 

The distance r^ shown in figure C3 is the distance along the curved path from 
the end of the turn when roll-out is started: 


^NOM 

ro = : Vg (Cl 6) 

^Limit 
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Figure C3.- Geometry of RNAV turn-anticipation and direct- to procedure. 


Direct-To Function 

When a direct-to message is issued to the aircraft directing it to turn 
immediately to C from its path AB, a calculation similar to that done for 
turn anticipation is performed- Equation (Cl 2) is executed in an iterative 
process to determine distance s in figure C3. In the iterative process a 
line is constructed from the point D* to the point C- The point D* is 
achieved by adding the roll-in distance r^ to the point D where the 
direct-to message occurs. Then equation (Cl 2) is executed using instead 

of A'F. This yields the first approximation for s. This value for s is then 
added to point D yielding a first approximation for E. The process is then 
repeated by constructing the line from E to C and continuing as before until 
the process has been performed five times. The purpose of this iteration is 
to obtain a good approximation for s. The function converges fairly rapidly, 
so five iterations were deemed adequate. 

A good estimate of s is needed because a direct-to RNAV constructs an 
imaginary waypoint a distance s from point D as shown in figure C3. This 
point is constructed so that the aircraft beginning its turn at D will fly 
to the path EC similar to the way it does toward BC in turn anticipation, 
that is, a smooth path having no overshoot. 
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METERING AND SPACING LOGICAL FLOW CHARTS 

Logical flow charts for the major M & S sequencing, scheduling, and 
schedule-maintenance control routines are presented in this appendix. Expla- 
nations of the speed-option eind DICE procedures are contained in the main body 
of the report. The logical flow charts which are included consist of the 
following: 

Figure 

Twenty-minute flight-plan message (FLTP20) Dl 

Five-minute flight-plan message (FLTP05) D2 

Master control routine for schedule determination 

and maintenance (TENTSBQ) D3 

Schedule time-of-arrival computation (SLTCOMP) D4 

Schedule slippage (SLIPSLT) D5 
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(a) Page 1 of FLTP20 subroutine. 


Figure Dl Twenty-minute flight-plan message (FLTP20) 
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(b) Page 2 of FLTP20 subroutine. 


Figure Dl Concluded 
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(a) Page 1 of FLTP05 subroutine. 

Figure D2.- Five-minute flight-plan message (FLTP05) 
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(b) Page 2 of FLTP05 subroutine. 


Figure D2.- Concluded 
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^ mum path ^ 
^qual to nomi- 
^\nal pathj^ 


OM time based on maximum 
path and speed control 
(XMDELAY) is the same as the 
OM time on the nominal path 
assuming slowest speeds from 
present position 


Compute OM time assuming 
maximum path and slowest 
speeds from present posi- 
tion (XMDELAY) 


Compute OM time 
assuming nominal 
path and nominal 
speeds (TETA) 


Compute an amount of de- 
lay time which will be 
reserved for future 
delays 


^^hould\^ 
naximum ET^ 
include hold 
^N^sibil 


The 

maximum 

OM time 

(XMDELAY) is 

the OM 

time 

(TETA) 

plus a 

maxi 

mum desi 

rable 

hold 

time 



(b) Page 2 of TENTSEQ subroutine. 


Figure D3.- Continued. 
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(g) Page 7 of TENTSEQ subroutine. 


Figure D3.- Concluded. 
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(a) Page 1 of SLTOOMP subroutine. 


Figure D4.- Schedule time-of-ar rival ccmputation (SLTCOMP) 
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(b) Page 2 of SLTCOMP subroutine. 


Figure D4.- Continued. 
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